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The large Ud theory is constructed for the metallic state of high-T c cuprates. It is based on the 
Emery three band model extended with the O x -O y hopping t pp in the Ud — > oo limit. The Ud — > oo 
mapping on the slave fermion theory is used. The time-dependent diagrammatic theory in terms of 
the Cu-0 hopping t p d starts from the locally gauge invariant nondegenerate unperturbed state with 
vanishing average occupation n d of the Cu state and builds a finite nd in higher orders. This theory 
is locally gauge invariant asymptotically, but replaces the d — p anticommutation of the fermions on 
the Cu and O sites by the commutation and is antisymmetrized a posteriori. Rather than t p d, the 
small parameter of the theory is Ud < 1/2. The lowest order of the U d = oo theory generates the 
Ud = single particle Dyson propagators of the hybridized pdp- and dpd-fermions which exhibit the 
covalent three band structure filled up to the appropriate chemical potential /i. The leading many 
body effect is band narrowing, different from that found in mean-field slave boson theories. It is 
accompanied by the broad incoherent background related to the dynamical quantum charge-transfer 
disorder associated to the d 10 -f->-d 9 Cu/02 intracell charge transfer fluctuations. The disorder effects 
fall well below the Fermi level and break the Luttinger sum rule for the conduction band. Those 
results and even the infinite order non-crossing approximation are insensitive to the omission of 
the d — p anticommutation rules. The contributions affected by the d — p commutation show up 
in single particle propagators beyond the third order. The effective local repulsion between the 
hybridized pdp propagators turns out to be a sizeable t pd /A d/1 , where is the difference between 
the energy of the Cu-site and fi. The a posteriori antisymmetrization of the theory removes the 
triplet repulsion between the pdp particles but keeps the singlet repulsion which favors at low energy 
the incommensurate SDW correlations. Such t pd /A d/1 repulsion is the metallic counterpart of the 
Ud = oo super exchange J p d between the dpd propagators. Resonant valence bonds appear thus 
as incoherent perturbative corrections here. The resulting modified slave fermion theory (MSFT) 
approximately obeys the local gauge invariance and conserves the local fermionic anticommutation 
rules, provided that rid is sufficiently small. The corresponding theoretical predictions compare 
favorably with ARPES, NQR, X-ray, neutron scattering, Raman, optic and superconductivity mea- 
surements, emphasizing the importance of oxygen degrees of freedom in the physics of high-T c 
cuprates. 

I. INTRODUCTION 

The long standing question in high-T c cuprates concerns the nature of interactions which are responsible for 
the superconductivity and the other unusual properties of those materials. The early high-energy spectroscopic 
measurements indicated that the Hubbard interaction Ud on the Cu site might be quite large. This opened the 
question of whether or not Ud alone can account for the basic physics of the high-T c cuprates. Such a question can be 
rephrased^ in terms of the structure of the effective interactions in high-T c cuprates, which include the concomitant 
strong electronic correlations. This fundamental question is discussed in some detail here from the theoretical point 
of view, with the results finally confronted to some salient experiments. 

The observed phase diagram of cuprates is characterized by a crossover (ignoring the small interplanar couplings) 
between the insulating long or short range AF phase at small hole doping < x < x cs to the disordered metal- 
lic/SDW/superconducting phase for x > x cs . Typical experimental values of x cs found from ARPE£&^ are of the 
order of a few percent. The local properties related to the measured ARPES spectra are the average single-particle 
occupations of the Cu and XiV sites rid and n p . The latter can be found in the x > x cs metallic phased from the 
electric field gradients measured^ by NQR, giving rid around 1/2 which increases on doping with holes. While an 
accurate evaluation of rid in the x < x cs regime is hindered by strong local magnetic fields in the AF phase and the 
narrow range of its stability, it is usually inferred^ that rid decreases slightly on doping. It is important to note that 
the NQR results in the well-developed metallic phase rely only on the measured local symmetry of the average charge 
distribution in the vicinity of the Cu and O nuclei. They are thus essentially model independent. 

The crossover is also clearly evidenced by transporter— optical^ and Rama n 14 ' 15 measurements. In addition 
to the low frequency conductivity, which presents an unconventional behavior for x > x cs , the strong optic edge, 
associated for x < x cs with the excitation through a gap, is smeared oufci^ for x > x cs into the transitions between 
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broad structures with finite threshold frequencies. Consistently, the Raman data exhibi t 16 i 17 a neat two magnon 
AF resonance for x < x cs that smoothens into a broad continuum for x > x cs . This observation directly reveals a 
deconstruction of the AF order upon doping. 

There are many other features of cuprates which also corroborate the proposed crossover picture. One such example 
is the incommensurate magnetic ordering. This ordering is associated with the magnetic superlattice Bragg spots 
at qsDW, which differ a 14 i 18 ~ — from qAF = G/2, where G — 2ir[l, l]/a is the reciprocal lattice wave vector (a is the 
Cu-Cu lattice constant). Furthermore, symmetry analysis shows that the incommensurate magnetic order may be 
accompanied by a % modulated charge transfer (CT's) within the Cu02 unit cell and the CDW among the cells, 
coupled in turn linearly to the incommensurate lattice deformations. The deformations give rise to the nuclear 
superlattic o 14 i 19 i 24 i 25 Bragg spots, with small wave vector q . Experimentally, qsDW and <fo are related through the 
simple Umklapp relation^ 1 - cfo + G — 2qsDW, irrespective of the value of x. When only one leg qsDW (and c/b) of the 
wave vector star is present the D4 symmetry is broken and the entangled 2 ^ magnetic, CT, CDW charge transfer and 
lattice deformations appea r 24 ' 25 in what is often called nematic stripesj 21 ' 23 For x < x cs , qsDW (and qo) lie along the 
diagonals of the CuC>2 zone (diagonal stripes) but rotate 1 - 5 - by 7r/4 for x > x cs to the positions along the main Cu02 
axes (collinear stripes). They show a commensurate LTO/LTT instability^ 7 . - — which gives rise to the new Bragg spot 
at 27r[l,0]/a that is coexisting with collinear stripes. According to ARPES data, the LTO/LTT lattice instability 
occurs for x rj 1/8 close to the dopin g 3 ' 30 ' 31 x v h > x c which puts the Fermi level on the van Hove singularity of the 
conduction band. This relates a prominent band feature, namely the vH singularity, to a phase transition observed 
for x > x cs - Remarkably, the LTO/LTT lattice instability suppresses^ 7 - the superconductivity in LBCO that would 
otherwise be close to optimal. In addition, the spin and charge disorders are observed in the glassy and metallic phases 
of cuprates by many experimental methods, including NMRj^Sr— transport) 11 ' 12 IR 9 -^ and Raman spectroscop y 14 i 35 

The simultaneous appearance of magnetic and lattice Bragg spots is a clear signature of the spin and charge 
coherence in the ground state of the system. The superconductivity, which occurs exclusivel y n i 12 i 15 for x > x cs , is 
itself a coherent state. Apparently, those coherent features are not only balanced among themselves but also compete 
with the spin and charge disorders. The main aim of the present work is to contribute to the understanding of these 
competitions. 

II. WEAK VERSUS STRONG COUPLING 

The theories of high-T c cuprates often start from the tight binding model, with the vacuum consisting of Cu + (d 10 ) 
and 2_ (p 6 ) states. The Cu 2+ (d 9 ) and 0~(p 5 ) site energies are denoted then by Ed, s px and e py with A p d = £ p —ed > 
in the hole language and A pp = e px — e py . The Cu 3+ (d 8 ) state is reached by spending the energy 2ed + Ud where 
Ud describes the bare interaction of two holes on the Cu site, which may be reduced to some extend by intra atomic 
correlations. In contrast to that the 0(p 4 ) configuration is usually associated with the energy 2e p i.e. U p is considered 
as relatively small. Referring 3 ^ to LDA results, such a model of the Cu and O sites was completed with the Cu-0 
hybridization t p d- This selection of the relevant single particle and interaction parameters is often called the Emery 
models The original Emery model was later extended^ 7 - by the direct 0-0 hopping t pp which describes the hole 
propagation rotated by ir/4 with respect to the Cu02 axes. The model is completed by choosing the total number 
of holes 1 + x per Cu02 unit cell, where x is the average number of doped carriers. The average single particle 
occupations of the Cu and O x , 9 sites nd and n p are then linked by the sum rule rid + 2n p = 1 + x. 

In order to explain the coherent features of cuprates, some early theoretical work a 1 ' 38 invoked a small Ud regime 
of the t pp = Emery model. This was later extended^ under the assumption Ud < A p d, to an effective intraband 
Ud reduced by metallic kinematics 4 ^ correlations. The x = Fermi surface touches the logarithmic vH singularities at 
T=0, irrespective of the ratio t p d/A p d- The equal sharing of charge nd — 1/2 at x = between Cu and two O's is 
obtained in the covalent limit s 38 ' 39 t p d S> A p d, whereas nd = 1 corresponds to the opposite ionic limit A p d 2> t p d- The 
x = Fermi surface is perfectly nested* 1 - ' 3 ^ 41 Finite Ud, and especially its Umklapp component^ enhances therefore 
quite strongly the commensurate qAF = G/2 SDW fluctuations. The latter may give rise 4 ^ to the unconventional 
behavior of the conductivity. The q = O x /O y CT fluctuations within the Cu02 unit cell are also enhanced, whether 
coupled 2 ^ or not^^ to the G/2 SDW. The q = Cu/0 2 and O x /O y CT's make the Raman activeii quadrupole 
moment of the Cu02 unit cell vary, while conserving^ 3 - its total charge. Therefore, in contrast to the optically active 
q — > CD W 42 ' 44 ' 45 or Raman active Ai g Cu-0 CT^ 4 - the q — > Bi g or B2 9 intracell charge (or current) fluctuations 
are not frustrated by the long-range Coulomb forces4^ When coupled linearly to the acoustic modes, the O x /O y CT 
may induce the acoustic lattice instability 2 ^ at q = 0. The q = O^/Oy CT also couples quadratically^ 3 - to the 
very slow tilting modes in lanthanum cuprates which result a 43 i 46 in the LTO/LTT instability of LBCO. With finite 
doping x the SDW instability move a 47 i 48 to incommensurate values of qsDW and becomes weaker in this model. In 
particular the effect of the G-Umklapp interaction Ud in the build up of AF correlations is diminished in this way. 
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Usually, the Umklapp in question is removed by hand from the theory when qsDW becomes incommensurate, i.e. the 
SDW commensurability pinning together with the resulting intrinsic striping and disorder is ignored. This results in a 
(too) smooth sliding of qsDW with x. The incommensurate O x /O y CT fluctuations show similar behavior, assuming 
that they are driven by two SDW's, which gives ris o 21 i 23 ' 26 to the observed relation % + G = IqsDW between the 
incommensurate lattice and magnetic wave vectors. Finally, in the presence of the attractive interactions between 
carriers, the weak-coupling prediction is L 38 ' 41 that the SDW instability is replaced by the superconductivity either by 
doping x or at x = by (chemical) pressure. Although the SDW behavior emerges more or less correctly, the main 
problem of this description is that it puts^ 3 . the commensurate LTO/LTT instability at x = rather than at sizeable 
hole doping x ps x v h ~ 1/8, where it is observe d 19 ' 27 ! 29 in LBCO. 

This could be remedied by including t pp . Although smaller than t p d on chemical grounds, it is quite relevant 
in the weak coupling theory. First, for t pp < (appropriate 3 - 7 -^ for high T c cuprates), t pp set o 37 ' 43 ' 49 the Fermi 
level of the x = half filled lowest band below the vH singularity, which means that the latter is reached upon 
a finite hole doping x v h > 0. The ARPES spectra of the hole-doped cuprates in the x > x cs metallic state can 
than be fi t 37 ' 49 by the Emery three-band structure. Those fits indicate^ 9 - that the bare parameters obey the relation 
A 2 d > 2tp d > Apd\t pp \ > At pp with \t pp \ large enough to account 3 - 7 - for the tt/4 rotation of the Fermi surface (Fermi 
arcs) with respect to its t pp = form. Concomitantly, the Cu occupation rid at x = is reduced below 1/2 in the 
covalent limit t p d > A p d- Furthermore, t pp break o 50 ' 51 the perfect nesting properties of the x = Fermi surface, i.e., 
using the Id language^ it plays the role of the imperfect nesting parameter. The elementary SDW particle-hole 
bubble develops then the peaks at incommensurate qsDW f° r x = x vH and small o;j 37 i 48 ' 50 When the small interaction 
Ud is introduced, the resulting G-Umklapp scattering of two particles is in discord with this value of the wave vector 
and the resulting value of qsDW = (<fo + G)/2 is, in general, incommensurate and weakly affected by a small Ud- On 
the other hand, the elementary O x /O y CT particle-hole bubble for x = x v h is2& logarithmically singular at u> = 0, 
q = for any value of t pp , corresponding to the J aim- Teller splitting- ^ 43 ' 46 of the vH singularities, and favoring the 
commensurate LTT instability. While this latter result agrees with observation a 19 i 27 i 29 in lanthanates, the problem 
is that in the weak coupling theory, strong magnetic correlations occur only for x w x v h- As mentioned above, 
this particular feature is at variance with observations where appreciable magnetic correlations coupled to the lattice 
(stripes) persist over a wide range of doping, from x > at least up to optimal dopings. While the weak coupling 
theory with finite t pp explains thus the metallic phase reasonably well, it fails to describe the Mott-AF phase for 

X <C X C s • 

Another important feature not encountered in conventional weak coupling theories is the intrinsic CT disorder 
d 10 +p 4 <->d 9 +p 5 (p n 's will be dropped from now on), as well as the corresponding d 9 spin disorder on Cu-sites. 
Already in the earliest theoretical works with Ud S> A pd , it was pointed out 5 - 2 - that the static d 10 <->d 9 disorder is an 
essential feature of the t p d = limit when Ed falls within the dispersive band, A p d < 4\t pp \ in the present language. 
A finite t p d is expected to render the d 10 <->d 9 disorder dynamic. Indeed, a broad dynamic background appeared 
in the early 5 - 3 - slave particle NCA calculation and also in the non-magnetic^ and magneti o 55 ! 56 DMFT calculations 
with Ud > A pd - However the relation between the broad background and the dynamic d 10 <->d 9 disorder had not 
been established until recently^ It is noteworthy that the strong ionic electron-phono n 43 i 52 coupling may make the 
d 10 <H-d 9 disorder slow again by a polaronic reductio n 33 ! 52 of single particle hopping and so account for the observed 3 - 4 - 
quasi-static intrinsic charg o 33 i 52 and spi n 42 i 57 disorder of cuprates, reflected in the unusual behavior of transpor t 11 ' 12 ! 42 
coefficients and NMR relaxation ratesi 33 ' 34 

All this motivates us to investigate carefully the Ud 3> A pd limit of the £ pj ,-extended Emery model, omitting 
at present the electron-phonon coupling (although it is possibly strong^ 3 -). In this endeavour, we are led to some 
extent by the translationally invariant Ud = oo mean field slave boson (MFSB) theory; 5 - 8 - which predicts for optimal 
dopings^ 9 - the band picture with weak renormalization of the band parameters t p d and A p d (t pp is unaffecte d 49 ' 5 1 1 59 ) 
in the regime A 2 d > 2t pd > A pt ;|ipp| > 4f 2 p , identified above on neglecting their renormalizations and keeping 
them fixed for a given parent material. At x — the t pp = MFSB describes 5 ^ the change of an insulator into a 
correlated metal through the Brinkman-Rice (BR) phase transition between the n d ^ FSB = 1 and < < 1 

states for (A p d/t p d)BR ~ 4.7 i 58 ' 59 This transition is conserved^ 9 - for small t pp but shifted linearly in t pp to higher 
values of {A pd /t pd )BR- In the MFSB language, the cuprates with A 2 d > 2t pd > A pd \t pp \ > 4t pp fall below 5 -^ the 
BR transition at x = (i.e. A pd /t pd < (A pd /t pd )BR)- For x finite, the MFSB smoothes out the BR transition 
in n^f FSB (x). Close below the BR transition a few percent doping of the x — state can then easily produce 5 - a 
sizeable decrease of n d AFSB from the x — value n d AFSB « 1, i.e. dn^ IFSB / dx < in the BR regime. Further 
below the BR transition, the renormalizations become weak and the weakly renormalized metallic dn d IFSB jdx > 
regime is obtainedj 5 - dn d IFSB /dx = conveniently defines the position {A p d/t p d) C s < {A p d/t p d)BR of the crossover 
and gives the corresponding doping value x^ F B (A p s d ,t p s d ,t pp ). Typically^ 5 - n^f F (x^ FSB ) is close to 3/4 already 
for t pp = 0. The x = value n d AFSB is, however, somewhat overestimated in the BR regime keeping in mind that the 
BR transition should itself be smoothed 4 ^ in a theory better than MFSB. In particular, the MFSB does not contain 
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the d 10 -H-d 9 disorder and the magnetic correlations. The latter give additional stability 6 ^ to the insulating phase 
and open the possibility of crossing over from the AF-insulating regime to the metallic regime. Moreover, the local 
gauge invariance, although satisfied on average in the MFSB,— is irremediably broken . 49 ' 62 Altogether, this leads us 
to search beyond the MFSB. 

Extension of the MFSB is usually attempted from the x = 0, Mott-AF side, starting from the unperturbed 
n^ ' = 1 Neel ground state. That state is widely used to approach the propagation of the first additional hole^ 3 - 
or electron^ in the x=0 AF phase of cuprates. The doped hole is place d 63 ' 65 on the upper oxygen level, assuming 
that rid, associated with the lower copper level, is close to unity. Originally, the t — J model with the large Ud 
superexchange J p d = 4t 4 d /A 3 d was so obtaine d 59 ' 63 ' 66 ' 67 for Ud > A p d 3> t p d > and t pp — 0. In the opposite limit 
Ud > tpd 3> ^pd > appreciable t p d hybridization within the Cu02 unit cells leads for x — to a reduction of nd, 
typically to « 1/2. It is then preferable to pu t 26 ' 68 the additional hole on the lower (rather than upper) covalent 
level, let its average charge in the unit cell be shared equally between the Cu and O states, worry about the single 
particle hybridization among the unit cells, ensure that holes on Cu sites avoid each other in the temporal dimension 
and allow for the spin polarization on Cu and O sites, as emphasized in the preliminary report 2 ^ of the present work. 
This means that the t — J model is not a suitable starting point for t p d 3> A p d, i-e. that the covalent phase, as 
described here with dynamic spin and charge disorder, including weak magnetic correlations, prevails then even at 
x = 0. 

The intermediary regime A 2 d > 2t^ d > A p d\t pp \ > 4t pp requires additional care. This regime is apparently at 
the brink of instability of the coherent x = Mott-AF phase. The crossover is then expected^ to occur already 
for small x s=s x cs . Nevertheless, the t — J approach is often rigidly extended to all x > of interest, taking that 
additional holes go to upper covalent levels upon tacit assumption that x < x cs in cuprates. Adding phenomenological 
next-to-next-Cu-neighbor effective hoppings t' ', t" . . . and, sometimes j* 9 - next-to-next-Cu-neighbor super exchanges J' , 
J" is obviously insufficient in this respect, essentially because intra- and inter-cell Cu-0 covalenc o 26 ' 70 of the lower 
level and the associated temporal incoherencies^ are omitted. The restriction to rigid t, t' — J, J' models can in 
principle be relaxed and the x < x cs theory based on the n d = 1 Neel unperturbed ground state of the Ud = oo 
Emery model extended to dopings x > x cs but this requires high order calculations. We are thus tempted to restrict 
the rigid t — J approach to the range x < x cs , including 6 ^ the required 7 - - covalent corrections in that limit, while 
here we describe the doping range x > x cs in terms of a renormalized three-band theory, determining x cs from the 
upper side. Such an approach bears some resemblance with nonmagnetic MFSB ) 49 ' 58 NCA ; 53 ' 59 DMFT— ~— and 
with LDA+U 7 -^ calculations for cuprates. It contains dynamic charge and magnetic fluctuations on Cu and O sites, 
including perturbatively the resonating valence bonding (RVB) related to the generalized superexchange J p d- 



III. U d = oo PERTURBATION THEORY IN TERMS OF t pd 

The T = diagrammatic expansion in terms of t p d is used for this purpose through the slave particle mapping, 
which is asymptotically locally gauge invariant. The time structure of the perturbation theory plays an essential 
role because Ud = oo is replaced by time delay s 26 ' 39 ' 40 ' 61 , induced in the motion of holes across Cu sites by the 
propagation of the intermittently added single particle. The corresponding n? = 1 slave particle theory is quite 
intricate, especially 6 ^ for t pp ^ 0, as also indicated by recent small cluster calculations^ Here we start therefore 
immediately with the metallic phase for appreciable doping. Partial infinite sums in terms of t p d are selected with the 
regime 2t 2 d > A p d|ipp| > 4t pp in mind. t p d is not the small parameter of those expansions but rather, on noting that 
the effects of large Ud are reduced if nd is small, the small parameter is rid < 1/2- The effective repulsion, which is 
proportional to t pd (as is the attractive superexchange in the opposite limit), turns out to be small to comparable to 
t p d, provided that nd < 1/2. The metallic nd ~ 1/2 regime is thus reached in reasonably low order of our diagrammatic 
Dyson perturbation theory, which gives a practical value to these summations. It is shown from the x > x cs side that 
x cs (A p d, t p d, t pp ) bears then some resemblance to x^ F (A p d,t p d,t pp ). We also show how the d 10 -O-d 9 disorder and 
the SDW correlations enter the single particle propagation, while the singular properties of the coherent O x /O y and 
SDW correlations, associated with imperfect nesting, already briefly described elsewhere^ will be further discussed 
within the present x > x cs approach in an upcoming publication^ 

The Ud = oo theory with auxiliary (slave) particles 7 - 2 ^— is well known and will be discussed only briefly here. The 
d 10 state on Cu at the position R is denoted by /g|6) and the d 9 state with spin a by ^ |6), where |0) is the auxiliary 
vacuum on Cu. In the so spanned three-state space (ds state at led + Ud omitted) , the number operators of the slave 
particles satisfy = n^ + J2 a n 1a = 1- The physical fermion c°? projected on the d 9 , d 10 subspace is written as 
c°J — > b a ^fg. The corresponding number operators satisfy = n^, usually called the Luttinger sum rule (LSR). 
b a J and ft can be taken respectively as fermions and bosons ("slave boson theory", SBT) or as bosons and spinless 
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fermions ("slave fermion theory", SFT) in order to satisfy the anticommutation rules on and among Cu sites projected 
on the dg, dio subspace. The states on oxygens are associated with physical p-fermions. The U d = oo Hamiltonian 
written in terms of the auxiliary particles is locally gauge invariant, and commutes with Qg, i.e. = 1 is a physical 
constant of motion. 

Therefore we start the time-dependent perturbation theory, in terms of Hi(t p d), from the unperturbed, t p d = 
0, Qfi = 1, paramagnetic, translationally invariant slave particle ground state associated with the unperturbed 
Hamiltonian H QX - XN = H Qd (e d ) + XY,r(Qr ~ 1) + H Qp (e p ,t pp ) 

\G x) = \G b ) ® \G f ) ® |Gg) , (1) 

where 

i G o>=n4i 6 ) (2) 

R 

is the = 1 state. The corresponding energy A is the site energy of the /-particle appearing, as usual, by adding 
X(Q £ — 1) into the unperturbed Hamiltonian Hq. In the present perturbation theory A serves as a parameter which 
checks the local gauge invariance, rather than as the Lagrange multiplier familiar from the MFSB theories. 

\Gq) is the state with no ^-particles, i.e. = rv d = 0. Thus n d ^ = is the outset of our expansion in terms of n d 
small. Since \Gq) is the no-particle state it is nondegenerate, irrespective of the Pauli symmetry of the b^-particles. 

The /-particles in rSp = 1 Eqs. (|H2p can also be chosen as bosons or spinless fermions. The advantage to choose 
/'s as bosons is that then 6 CT 's can be taken as fermions indistinguishable from p- fermions. In the slave particle 
representation this permits one to satisfy the anticommutation rules between the physical c-fermions on the Cu 
sites and the p-fermions on the O-sites. The disadvantage is that the state ([2]) with bosons is highly degenerate with 
respect to multiple boson occupations of Cu sites. As is well known, degeneracy of the unperturbed ground state causes 
problems in time-dependent perturbation theory. This difficulty is obviously eliminated on choosing fa as spinless 
fermions and fo^'s as bosons, since the state of Eq. ^ is then nondegenerate. However, in that case, we are dealing 
with three kinds of distinguishable particles, i.e. the anticommutations between the Cu sites and O-sites are replaced 
by commutations. The corresponding time-dependent perturbation theory (SFT) must be therefore antisymetrized a 
posteriori (it will be then named here "modified" SFT, MSFT) and this is the route chosen henceforth. 

As the state of Eq. @ is nondegenerate for spinless fermions it can be simply expressed in terms of Fourier 
transforms ft of the local operators ft. Up to an unimportant phase factor we have 

iG£>=n4i6>=n$ 5 >' (3) 

R k 

where the product over k extends over the whole Cu02 Brillouin zone, which corresponds to N Cu02 unit cells. In 
other words, the Mott state of spinless fermions is equivalent to the full (dispersionless) band of those particles. 
Let us finally mention that \Gq) is the usual (nondegenerate) Hartree-Fock (HF) state of the p-fermions in the 

cosine band associated with e p and t pp . This band contains 2n p ^ = 1 + x fermions, associated with the chemical 
potential For x small this band is nearly quarter filled. It is usually folded artificially into the Cu02 Brillouin 

zone in two 1,1 oxygen bands, anticipating the effect of Hi(t pd ), which is expected to generate three separate bands 
and make the lowest one nearly half filled. 

Once the unperturbed ground state is expresse d i 53 in the full momentum representation so should the slave particle 
Hamiltonian H\ = H \ — AA^ + Hi. We have, in terms of ft, bt , p^ (i = I, I) and their hermitean conjugates, 

i,k,a k,a k 



(4) 



t { M) 



'|sin^|±|sin^ 
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assuming the D4 symmetry. Here, i p j (fc) describes the fact that by annihilating the f$ spinless fermion and by creating 

the bt^^ boson, one annihilates the fermion in either of two i = 1,1 bands = e p ± 4|i pp | sin (k x /2) sin (k y /2). 

A is introduced in Eq. ^ as the test parameter (rather than as the Lagrange multiplier) since the physical = 1 
result, independent of A, must be ultimately achieved. As will be seen below the present time-dependent perturbation 
theory will prove independent of A in each order. 

The perturbation theory can then be carried out in terms of Hi on the top of the nondegenerate state of Eqs. (|f I2I3[) . 
the time orderings and the normal orderings being well defined together with the Pauli symmetry of the relevant b-, 
/-, and p-particles. Since both the Hamiltonian H\ and the unperturbed ground state are locally gauge invariant, 
translationally invariant on the Cu02 lattice, and symmetric under time reversal, the SFT will either generate the 
exact ground state with the same symmetries, or break them in a controlled way. It is thus left to the SFT to keep 
Qj^ = 1, generate the LSR rid = ?^t>, obey the anticommutation rules on and among the Cu sites, and to satisfy the 
charge conservation rule rid + 2?i p = 1 + x. Such a multiband SFT does not suffer from problems related to the 
breakdown of local gauge invariance, encountered^ in the single band models. The MSFT is eventually constructed 
only to take care of the Cu-0 anticommutation rules, as well as possible. 

We wish to emphasize the intimate relation between local gauge invariance and causality in the time-dependent T = 
perturbation theory. The latter has therefore advantages over the finite- T Matsubara theory often encountere d 53 i 67 
in the slave particle context. The Matsubara theory uses the canonical ensemble in the full slave particle space and 
therefore treats the = 1 and the 7^ 1 states on equal footing, provided that they are degenerate in energy. It 
is therefore more difficult to control the Q ^ = 1 local gauge invariance in the finite T Matsubara theory than in the 
time-dependent T — approach. 

The elementary bricks which build the time-dependent perturbation theory according to Wick's theorem are the 
free-particle propagators. Defining, as usual, B\(k, t) = —i(Tb^bt(t)), we find that the free propagator of the 6-particle 
is dispersionless, 
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u — £d — A + ir\ 



*v = - — ■ (5) 



Through +ir] it describes the intermittent creation of the ^-particle, while its annihilation is impossible in the no- 
bosons state of Eq. (fTJ). In contrast, the spinless fermions can only be annihilated, i.e. F { x 0) is obtained from Eq. ([5]) 
by replacing Ed + A by A and +ir/ by —ir\. The free propagators of 1 + x p-particles contain both +ir] and — irj 
components Gp > {k,u) and Gp < (k, u>) according to their Fermi distribution in the HF state associated with the 

i = 1,1 bands. The corresponding chemical potential is hereafter denoted by With x < 2, only the states in the 
/-band are occupied. 



Figure 1: (Color online) Propagator D^ ' of the spinless fermion (green)-boson (blue) pair carrying wave vector k and energy 
w; arrows denote that the 6-particle (blue) can only advance and the /-particle (green) recede in time. 

The d-particle propagator is defined by D^(t) = —(i/N)(TJ2 q f\bk+qfq{t)b\ +q {t)) . shown diagramaticaUy in 

Fig. [U is thus also dispersionless and obtained from Eq. ([S]) by replacing Ed + A by Ed- = (oj — Ed + ir))^ 1 thus 

firstly reproduces n^ ' = and secondly is independent of A, as any physical quantity should be. In the next step we 
define D^S associated with the r-th order time-dependent perturbation theory r > 0. According to the definition of 
a particle-hole pair, D^'(lj) is given by the (generalized) Bethe-Salpeter equation 

D f H = ^ H + sf- 1 ' ( W )r<r> HE<f H . (6) 

Here sl r 1 ' is the quantity irreducible with respect to cutting the p-lines and rl' ' (cu) is the renormalized four-leg 
vertex given iteratively by the Dyson equation 

rfH = rf ( w ) + rj?> MEjr%)r< r >( w ) . (7) 
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in terms of the bare four-leg vertex ri ^^) characterized by fJ,^°\ shown in Fig. [2 

4° } M = [t^{k)] 2 Gf<{k,u) + £[^(fc)] 2 G«>(fc, W ) . (8) 



i=l,l 

Eq. (JHJ) generalizes the four- leg vertex used previously^ for t pp = 0. 



Figure 2: (Color online) Four-leg vertex (uj); triangular vertices are t p ^(k) and the red lines are the free propagators 
Gp (k, uj) combined according to Eq. {5). 



Actually, t'^rV in Eq. ft) can be interpreted as an appropriately symmetrized generalization to the Emery model 
of the wide-band propagator on the Anderson lattice. According to Eq. ([5]) the intermittent p-particle t~ d T^ is 
prepared in two i = Z, / p-bands, instead of one. t~ d T^ is thus the canonical pdp propagator which appears naturally 

in the perturbation theory for the Emery model. For r > 0, this propagation involves both i = 1,1 bands and the 
d-state, similarly to Eq. (JSJ) for an intermittently created d-particle. Eq. ([7]) can thus be interpreted as the Dyson 
equation for the t~ d T^ single particle pdp propagator with the Dyson self energy Il' r ' = t^ d ^S ^ irreducible with 
respect to t~ d T^ -lines. According to Eq. (0| for Hi, the lowest order si ' is simply equal to of Fig. [T] In 
other words is not only the elementary d-particle propagator but also the essential component of the lowest 
order "local" irreducible self energy II' 1 ) = t^ d D^ in Eq. ([7]) for t~ d T^ . The r = 1 procedure thus separates out 
the ^-independent free d-propagator in the leading pdp-particle self energy II^ 1 ) on associating in Eq. ([5]) the 
fc-dependence of triangular vertices of Fig. [2] with the fc-dependent free p-propagator t~ d T^ ' . 

On the other hand D j (u>) describes the creation/annihilation of the intermittent e?-particle on the Cu sites and 
its subsequent dpd propagation. The factors Y^S ^(uj) in Eq. © are the same as those involved in rl r ' ) . Due to this 
sl 7 1 ' can be taken as the effective "free" propagator on the Cu-site, which allows one to perform the resummation 



of the Bethe-Salpeter Eq. (JB|) into the Dyson form, in accordance with the idea that DS is the propagator of the 
single dpd particle. Nevertheless, sl r ^ is not in general a fermion propagator and, consequently, neither is D*j?\ 
The reason is that the SFT removes the d 8 state completely, while the original large Ud theory treats it as the empty 

(r) (r) 

state (upper Hubbard band) at large energy Ed + Ud- Thus, in order to satisfy the equality n d ! = n d _ associated 



should allow for the additional spectral density at the energy Ed + Ud on the t < side. This adds a term 



with the full fermion anticommutation rule on the Cu site and that obtained from D V (t) in the t — > 0± limits, one 



(r-l)> _ 



a 



(r-l) 



uj - (£ d + U d ) + in 



(9) 



to E^ -1 ), which takes explicitly into account the fact that the d 8 state is empty. Although permanent (average) 
occupation of the d 8 state is forbidden in the Ud — > oo limit, the full fermion nature of e?-particles requires visits of 
the d 8 state on the t < side. The point is that when the expression (|9|) is integrated over uj in the Fourier transform 
which determines the overall spectral density associated with the t < component of E^" 1 ) it gives a contribution 
a( r ~ 1 ) which has to be retained in spite of Ud — > oo. In principle, a lj '~ 1 ^ may be determined from the anticommutation 
requirement n d r }_ = n d _ of the extended Eq. © and is expected to be small for n d ) small. In practice, following 
the spirit of the slave particle theories, we associate the physical average occupation n d with n£: of Eq. ©, not 
worrying about o'( r ~ 1 ) > at all. 

The last step is to determine the chemical potential // r ) of the p-fermions in order to satisfy the conservation of 
the total charge 1 + x. In the SFT this can be done by requiring + 2n p r ^ = 1 +x, with the number of bosons n^ 

found below. Alternatively, one can require n d + 2n p r ^ = 1 + x on keeping in mind that the LSR nb = nd, though 
approximate in low orders, holds asymptotically in the exact SFT. 
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IV. LEADING CORRECTION TO SINGLE PARTICLE p- AND d-PROPAGATORS 



Here we choose the nj + 2n p = 1 + x prescription because, as we shall see now, Eqs. ([6][8]) at r = 1 then become 
equivalent to the HF equations for the free propagators in the tpd hybridized Ud — model. The deep reason for this 
far-reaching result is that the Cu site is initially empty, Eqs. ([T][3]), and at the order r = 1 the intermittent particle 
does not probe two-particle effects which involve the Cu site. 

More formally, Dt involves T,^ = = (w — Ed + i^) -1 , which is associated with unspecified commutation 

properties through the absence of the — ir\ component. The single particle problem of anticrossing^ 9 - between the Sd 
level and the "two" oxygen i = 1,1 bands is then solved exactly by Dg and rl 1 ^ of Eqs. ([B]|H]). Both t~ d T^ and 

in the Dyson form exhibit coherent poles belonging to three bands (branches of poles) oj^ denoted respectively by 
j = L,I, U. The poles wl in the HF t~ d T^ are associated with the residuals (spectral weights) z^\uj^) which can 
be expressed entirely in terms of the three . For example, for the lowest band L 



^\^ ] ) = — 1 fc m p " k m pk , (10) 

k y k ' .2/. II) , 1U) (L)s ' V ) 

and similarly for other two bands. The chemical potential fj,^ 1 ' of the p-fermions is next defined as the energy which 

— 2 (1) 

separates the poles of the p-propagator t pd Tj. in the upper and lower w-plane. Since the same poles appear in 

dQ\ this step accounts for them too, with no reference to the Pauli symmetry of the i>t/ pairs. In contrast to fi^°\ 
which defines the average number of p-fermions on the O-sites and allows for their fluctuations among those sites 
with the total number of p-fermions fixed, fj,^ allows also for fluctuations of the total number of p-fermions by their 
conversion into the f pairs. fiP-' {A p d, t p d, t pp , x) is determined through the approximate charge conservation rule 

n d l) + 2n p 1] = 1 + x, bearing in mind that rip and are defined by Eqs. (|6I7|) as functions of the band parameters 
and fj,^\ The whole r = 1 procedure described above thus amounts to the redistribution of the spectral weights and 
the Fermi occupation factors jy (with accompanying ztijy's) from two oxygen bands i = 1,1 and the empty d-state 

into the three hybridized bands uii j = L,I,U. In other words, the unperturbed ground state of Eq. ([T]) evolves 

through the prescription n d 1] + 2n { p ] = 1 + x into the HF state of the coherently hybridized noninteracting (Ud = 0) 
pd particles with a shift p« - in the chemical potential from the upper to the lower hybridized hole states. 
The shift is large when single particle anticrossing is important. However it is immediately evident that local gauge 
invariance is not obeyed for r = 1, because the double occupation of the Cu-site is allowed in the pd hybridized HF 
state. This will be corrected in higher orders. 

Concerning the anticommutation rule on the Cu site, we can take the t — > limit of (t) in the Dyson form, 

to find that n d _ = n d ) = n^, i.e. = in Eq. ©. The anticommutation rule on the Cu site is thus satisfied. 
However it is immediately evident that local gauge invariance is not obeyed for r = 1, because the double occupation 
of the Cu-site is allowed in the pd hybridized HF state. Later we shall return to this point more formally. 

Although the above argument, which shows that the Ud = oo SFT generates the Ud = HF result in the lowest 
order, does not require explicitly that ti^ 1 ' is small, the SFT, which starts with nf^ = 0, will converge quickly to 
satisfy local gauge invariance only when this condition is met. Let us therefore mention briefly the values of the 
single particle parameters A p d, t p d, and t pp which make n d (A p d, t p d, t pp , x) small for a given 1 + x. These conditions 
can be taken over directly from the three band HF theory^ which determines rS^ from the partial derivative of the 
HF energy of Eq. (Q} with respect to Ed, thus circumventing the clumsy calculation via the spectral density of -D^ 
complementary to that of Eq. (fTT))) . 

The simplest situatioi>L i 38 i 39 corresponds to t p d S> A p d and t pp = 0, where one immediately finds n d =1/2 taking 
formally x = (having the metallic regime with small x > x cs in mind). Note, in this respect, that /i^ coincides 
for x = with the vH energy lu v h in the lowest L-band and that the latter is almost independent of Ed in the 
limit considered. This results in the equal sharing of the charge between one Cu and two O's which, indeed, was 
traditionally obtained^ i 38 i 39 in this way. 

Finite t pp can be easily included in this scheme perturbatively for \t pp \ <C t p d, A p d- In contrast to t p d, t pp shifts the 
Fermi level fjS 1 ^ at x = from the vH singularity at oj v h- A finite doping x = x v h is therefore required to reach the 
vH singularity. x v h was founcl^ 9 - to be equal (t pp < 0) to —32t pp /n 2 A p d in the limit A p d > t p d and the corresponding 
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HF energy was determined analytically. This can be readily extended to the t p d > A p d limit with —t pp A p d <C t 2 pd 
when x v h ~ —32t pp A p d/Tr 2 tp d . The large numerical factor 32, multiplying t pp in this equation, is due to four t pp 
bonds per two t p d bonds in the CuC>2 unit cell. This compensates, to some extent, the chemical inequality t p d > \t pp \- 
The analytic calculation^ of the HF energy follows the same lines. A small \t pp \ is thus found to reduce the value of 

n^}\ In particular, for t p d ^> A p d the x = value of rvP falls below 1/2. 

These results can be further extended to the physical regime^ with sizeable t pp , satisfying 2t 2 d w —t pp A p d. This 

regime also sets (jp) in the L-band, while, importantly, the /-band is flat at e p . t pp is quite efficient in reducing 
below 1/2 at small x, due to a large numerical factor carried by t pp in this regime as well. Although the roles of t p d 
and t pp are comparable in the dispersion, a sizeable t pp requires, in contrast to t pd , a sizeable doping x — x v h to 

reach the vH singularity. In other words /Lt^ 1 ) for x = falls well below u) v jj (see e.g. Fig. [5k), i-e. n d is reduced 
appreciably with respect to 1/2. 

The useful overall rule of thumb is that n d < 1/2 as long as /Z 1 ) falls below the vH singularity at uj v h (independent 
of tpp), i.e. as long as x < x v h- In other words, the condition fii)? rj uj v h can be used to define the lowest approximation 
x^J (Apd, tpd, tpp) to the crossover function x cs {A p d, t pd , t pp ) although dnp jdx > all over the parameter space of 
interest. 



V. APPROXIMATE LOCAL GAUGE INVARIANCE AND THE LSR 



Here we turn first to the properties of the propagators B^\k,ui) and F^\k, ui), required to asses the accuracy 
of the above HF result and to construct the next order r = 2 iteration of Eqs. (|6l7l9p . The lowest order connected 

diagrams for AB^ (k, uj) and AF^ ] (k, w) are shown in Figs. 13141 As usual, the disconnected diagrams, which describe 
here the incoherent d 10 O d 9 vacuum fluctuations unaffected by the intermittent particles, are not shown. The arrows 
of time, associated with ±irj factors in elementary propagators B^\k,io) and F^\k, w) of Eq. ^ (independent of 
k) are depicted in order to emphasize that it is important to account for the full temporal structure of the theory. 
The bubbles which appear in Figs. [3] and [4] are the lowest order irreducible Dyson self energies for B^ (k,u>) and 
F^\k, lu). They both contain rl ^^) of Eq. ([8|) and involve summation over the occupied states in the /-band as 
indicated in Figs. I3|4I by the p-propagator going (only) backwards in time. In Fig. [3] the external frequency enters 
the bubble from right to left, while in Fig. |4]it goes from left to right. In addition, the spin factor 2 multiplies the 
/-bubble, unlike the 6-bubble, but we remember that there are two 6-bosons for each k. This symmetry, exemplified 
by Figs. 13141 will be referred to here as the / f-> b symmetry. 



Figure 3: (Color online) Lowest order renormalization for B. The arrows of time are shown. The energy and momentum enter 
the bubble from right to left. 




Figure 4: (Color online) bo/ symmetric lowest order renormalization for F. 

It means that the two diagrams of Figs. 13141 taken together keep the average value of close to unity. However, 
as easily seen, AB^(k,uj) and Ai^fc,^) are singular when the anticrossings^ between the e d level and bands 
i = 1,1 occurs for the states occupied in the /-band. This leads to a large difference between n]p of Eq. ((6]), which 
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takes the single particle anticrossing into account exactly, and the number n b of 6-bosons, which corresponds to the 
Schrodinger perturbation theory for the occupied single-particle states. The resulting convergence of the SFT is poor. 

It is therefore essential to re-sum the perturbation theory by associating the receding "oxygen" line in Figs, til II 
with hybridized rather than with For example, the 6-bubble is ^-independent ("local"), and given by 



JV - — UJ+Ujl '+X + ZIT] 

where are the residues of the propagators t~^T^ given by Eq. (fTU)) . It is also assumed for simplicity that t p d is 

large enough to have^ only the states in the lowest L-band occupied. Eq. (fTTj) is obtained by taking into account 
the Pauli nature of the involved particles, the energy conservation in each triangular vertex and by noting that both 
propagators in the &-bubble are running backwards in time. The set of occupied poles in Eq. pip lies in the upper 

w-half-plane, which makes Ani , associated with AB^' of Fig. [3J finite for 1 + x ^ 0, independently of A. 

We are now in a position to discuss the LSR = rid to the first order. Although n} is given by the HF theory and 
n^ 1 ' is evaluated in the Appendix it is more instructive to carry their comparison as follows. AB^ given by Fig. [3J 
with pdp hybridized rl 1 ^ on the "oxygen" line gives An£ . can be conveniently determined from the second term 
vj(o)pH-) j](o) £ ftx e Bethe-Salpeter Eq. ^ before its resummation into the Dyson form. On noting further that the 
product of the squared pole and a single pole of such Eg occurs also in AB^ , one finds easily that 



At-4 1 ' = n d 1] . (12) 

Although Eq. (fT2")l departs from the LSR rib = n<i, the difference between An£ and rvj/ is apparently small if 

is small, as shown in detail in the Appendix. The higher order corrections to n£ n^ 1 -* can then restore the LSR 

quickly, as the SFT converges from Qg m 1 towards the local gauge invariance Qb = 1. 

Indeed, a conclusion similar to m rip holds for the average value of Qj^. A result analogous to B^\oj) is 

obtained for F^\u>) by / •<-> b symmetry. (jr^'(cj) has a similar structure as complex conjug ate of-$\u>) of Eq. (fTTjl 
except that a spin factor 2 multiplies the sum over the occupied states in the L-band. The exact equation, which is 
analogous to Eq. (fT2"]l . 



An { p = -An { b 1] , (13) 

is therefore transformed into the approximate relation w 1 for n b w Anj, 1 ' = small, as further discussed in 

Appendix. m 1 is then subject to higher order corrections which asymptotically enforce the local gauge invariance 
= 1. The discussion of Eqs. (|12I13[) uncovers thus the mechanism of achieving the local gauge invariance in the 
SFT. The point emphasized here is that physically significant results are obtained already in the low order SFT 
provided that the average Cu charge rid (but not necessarily the total charge 1 + x) is small. 



VI. BAND NARROWING AND d 10 «d 9 CHARGE-TRANSFER DISORDER 



Once B^\uj) and F^\uj) have been determined, they can be used to calculate E^ - (B^ ] * F^) in Eqs. ® and 

([7]), i.e. to advance the iteration one step further to find d! 2 -* and rj?-*. Using the relation B^\to) = Bq(w — \) which 

follows from Eqs. ([5]) and (jTTj) , and b f symmetrically F^ (uj) = Fq ^ (uj — A) , we find (after integration over uj — A) 

that EW is independent of A. It is further shown in the Appendix that, according to Eq. (|TT|) . B^ (uj) can be written 

in terms of one pole in the negative w-half-plane and a set of poles in the positive w-half-plane, B^ = B^^ + B^ K , 

and / hJ symmetrically for F^' = i* 1 ! + (^ ne superscripts < and > denote arrows of time). As usual (see 

Eq. (|A.2[) ) the relevant contributions to the convolution E' 1 ) ~ (B^ * -F^ ) come from the poles on the opposite 
sides of the w-axis, 
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E« = 3(B«*2f>) 



g(B«>*2f)< + fl«<*^>) 

= £(!)>+ s^X . (14) 

E^ 1 ' is "local" (dispersionless) in the direct space, i.e. it "occurs" on the Cu site. As shown in the Appendix, we 
obtain in this way (independently of A) 

E (i)> = v 6 / i j (15) 



w - e d + lit] 

E (D< 



tL 



N 2 



x 



V f { -. L] f { - L) A k',k-^ (16) 

Z-~< J k> J k" .. ~{ L ) ~(L) , _ \ ±V > 



rifr rvp , id, Ap p,, oj^p . , ti^Jp, in Eqs. (O an d (|T5|) can all be expressed in terms of the spectral weight zp^ of 

Eq. (fTU)) and the chemical potential fi^\ E.g., e^ 1 ' = £d + 3/3^ (w = + A) where Px(oJ = Sd + A) is negative 
according to Eq. Importantly, = E (0) = £> (0) for 1 + x = 0, i.e. the propagation of the intermittent particle 

in the empty band is free, subject only to the single particle hybridization. t pd in Eq. (|16p is exhibited in order to 

stress that the leading term in the t p d expansion of E^X is quartic in t p d , in contrast to that of E^ 1 ^, which is 
then quadratic in t p d because rip and are. This clearly demonstrates that I]' 1 ) of Eq. (| 14[) is not a fermion 

propagator, vp , np as functions of A p d, t p d, t pp and x, as well as the double sum in Eq. (|16[) . are evaluated in the 
Appendix in the N — > oo limit, where the sets of dense poles are treated as cuts in the w-plane, once the temporal 
decomposition of F^' and E^ is properly determined. In particular, the shifts of dense poles can be neglected 
in the N — > oo limit, = 0J^\ = w- L ' in Eq. (fT()|) . We emphasize that the present theory, where d = 2 
explicitly, is not an expansion in the number of dimensions or in large orbital or spin degeneracy but in terms of rid 
small. Even the N ~ > oo limit is unessential, used only for the explicit calculation of the coefficients rip , np, id, 
A^p, Cjf ] in Eqs. dUI) and {TBJ). 

Apparently, E^ 1 ^ of Eq. (| 15[) has some features of the mean-field slave boson theory (MFSBT), namely in Eq. ([7]) 

(2) 

for rr ; it gives the band narrowing and the renormalization of the CT gap, 

t 2 -> t (1)2 - t 2 h 4- lr7 (1) ^ (1) 



A pd -34 1) (o;=4 1) +A) 1 (17) 



since e p , as well as t pp , remain unaffected. For rip small, when the LSR relation rip ps rip holds well according to 
Eqs. (dU), and QW = nP + rip ps 1, this reduces to t 2 d — > t 2 d (l — rip /2). Such t 2 d renormalization is about half 
of that predicte d 49 ! 58 by the MFSBT. Concomitantly, and in contrast to the MFSBT, is somewhat increased 
with respect to A p d- It is further interesting to extrapolate Eq. (|15p to larger A p d when np(A p d, t p d, t pp , x) ps 1/2. 
nP ps nP ps 1 - n p then gives t 2 pd {l + n[ 1) /2)n^ 1) closer to t^(l - nf FSB ) of the MFSBT. The difference here 
is that MFSBT replaces nP in Eq. (fT7| by n d IFSB , which connects it more closely to a self consistent than to the 
iterative r = 2 diagrammatic theory. Such differences are small for rip' small and the behavior of Eq. (fl7|) with 
increasing is satisfactory. 

We proceed by taking into account E 1 - 1 ^ given by Eq. (fT6|) . In contrast to the coherent band narrowing associated 
with E*- 1 ^ of Eq. (fT5|) . the receding (—irj) continua of Eq. (TTS|) describes the dynamic d 10 -f-)-d 9 CT disorder. Indeed, 
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the convolution in frequency means that the transformation of /^-particle (d 9 state) in the 6^-particle (d 10 state) 
is causally correlated in time. In contrast to that, the two k, k' integrations in £W< are independent, rather than 
being a convolution in k, k' = k + q. The "local" d 10 O d 9 CT disorder associated with Eq. ([TrJf is thus the result 
of a causal temporal and spatial combination of the "local" slave disorders appearing in "local" and which 
appear for 1 + x 7^ 0. These correlations are related to the local gauge invariance "to quartic order in t p d" , which 
localizes the CT d 10 od 9 event in a single Cu02 cell (it should be however kept in mind that the SFT is not locally 
gauge invariant unless it is carried to the infinite order in t pc j). E*- 1 ^ is thus the first step of the perturbation process 
which ensures that the /-particle is annihilated/created simultaneously with the creation/annihilation of the 6-particle 
in the description of the "local" CuC>2 CT disorder. In the limit TV — » 00, ImE' 1 ^ is a step-like function, finite in 
the range 2/iW — Sd > uj > 2ljm — £d- This disorder is revealed by the physical single particle propagators and 

t~dTf of Eqs. © and 0, defined by the irreducible self energy E' •'. Those propagators thus describe in particular 
the dynamic d 10 f>d 9 CT disorder of the permanent particles as seen by the intermittent particle. 




r x m r 







: ,v, 






r x- 



Figure 5: (Color online) (a) Im£ pd 2 rl 2 ' spectrum for A p d/\t pp \ = 8/3, i p d/\tpp\ = 5/3, t pp = —0.3 eV. Experimental points 
correspond to x = 0.07 LSCO; (b) M < 2 > is taken in Sec. IXI Al to satisfy the Luttinger sum rule. The leading d 10 -H»d 9 disorder 
effects are included schematically by gray and red shadings; (b) Characteristic to = /i*- 2 ' Fermi arcs obtained for < x < x v h 
correspond to the 7r/4 rotation of the Fermi surface induced by t pp . 

Let us thus consider of Eq. ([7| in some more detail, keeping in mind that a parallel discussion exists for 
The result is particularly transparent in the approximation which takes into account that for rA small ImE' 1 ^ is 
small with respect to the energy distance between the affected pole at oj^ and the Fermi level. In this, "quasiparticle" 
limit ReE^ 1 ^ is unimportant. The poles and the associated residue are then given by the coherent band narrowing 
lu^I of the three bands u>^ \ shown in Fig. O which is given essentially by replacing E' ) by E^ 1 ^ in rl 1 ^. z^ 1 - are 

j , k k k L ,fc 

the corresponding residues, such as those of Eq. (jTU)) . ImE' 1 ^ is then the only contribution to E' 1 ) beyond the 
coherently renormalized hybridization. In the relevant frequency range this gives 
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Imt^ifV) « V ^ . (18) 

p k V(^-^S) 2 +(^+^ ims(i)< ) 2 

where 77c = 77 sign/// 2 ) — u/ 1 !) and ImX^ 1 ^ is taken at w = u; ~. 

The step-like continuum in ImS' 1 ^ of Eqs. (fTr?|) and (fT8")) falls in the range 2wj\/ — Ed < cj < — e^, where wjf 
denotes the bottom of the L-band^ i.e. the disorder is far from the Fermi level. E^ 1 ^ is a receding term and, as a 
consequence, ImE' 1 ^ carries the same sign as 77s- This contrasts with the -F\-spinless fermion where the appearance 
of the disorder is related to the double change in sign of the imaginary part of the self energy (see Appendix) and 
so to the annihilation of /-fermions. Note that the sign of ImD' 1 ' could not be determined from the fact that the 
step-like continuum falls below the Fermi level but rather, its receding {—irf) nature had to be demonstrated directly 
in Eq. (jA.lOjl from its position in the upper half-plane. This illustrates the importance of the causal time-dependent 
analysis, carried out here. According to this analysis the continuum in ImE' 1 ^ is to be interpreted as the lifetime 
effect of the coherently hybridized occupied states, 

27T , „„v, (1) 



Pj-^^Im^X^') (19) 



rather than as related to the average p — d CT. In the quasiparticle approximation the effect of ImS' 1 '* is thus a 
broadening of the affected Dirac functions into the Lorentzians normalized to unity, shown by the broad line in Fig. [5] 
One should however take into account that, rather than zero, ImS' 1 ^ of Eq. (jA.lOjl is finite around the L^-band, 
all over the Brillouin zone in the frequency range 2u>m — £d < w < 2 ( u' 1 - 1 — Ed, where the coherent hybridization is 
absent. The coherent and incoherent frequency ranges contribute essentially additively to average occupations of the 
O (and Cu) sites. The conservation of the quasi particle spectral density in Eq. (fT9|) means in particular that ImE' 1 ^ 
does not affect the contribution rip HF ^ of the extended states to n p = n p 2HF ^ + rip for a given chemical potential 
/i' 2 ). Furthermore, by using Eq. ^ for Dj? either in the Bethe-Salpeter or in the Dyson form we can easily convince 
ourselves from Eqs. © and (fT5)) that, in addition to 



n^^n p 2H ^+n p 2 "^ , (20) 



we have 



n 



~{l + \n?)n?nr F) +n^, (21) 



d 



where n p 2mc ^ and n d 2mc ^ describe the incoherent contributions of ImE' 1 ^. In the simplest approximation n d 2mc ^ 
is obtained by integrating ImE' 1 ^ of Eqs. (|16p and (| A. 10|) over cj, while np inc ^ is neglected. The first term in 
np is somewhat smaller than n d 2HF ^ calculated from the HF bands defined by the self energy ImE' 1 ^ divided by 
(1 + /2)rvp , while the second term increases n d 2 \ The chemical potential fx^ is finally determined from the sum 

(2) (2,) _____ _____ 

rule = 1 + x. Equation (|2"Tj) combined with Eq. (|17j) establishes so the connection between the r = 2 

expansion and the familiar HF theory. 

For n^ 1 ' sa n^ 1 ' sa 1 — small Eq. Q2ip can be used, on linearizing in terms of n d 2H ^ — to write 



i) 



(2H) + 2n< 



2H) 



P 



X eff 



+ I n « 2 _ [nf nc) + 2nf n ^) . (22) 



**It - ■ 2 

corresponds thus to the HF chemical potential for the effective doping x e ff put into the L^-band renormalized 

according to Eq. (jTTJ) . n d 2H ^ can be then found straightforwardly, using the conventional HF procedure described 
briefly in Sec. IIVI Notably, Eq. (|2"2"|) explicitly breaks the standard Luttinger sum rule (to be distinguished from the 

LSR rid = nb)- Apparently, the shift /i' 2 ) — fjP^ which corresponds to the shift n d — is small for r^p small. 
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Indeed, the large shift of the chemical potential /j,^ 1 ' — fi^ due to the anticrossing of the d-state and two p-bands, 
obtained under assumption that n]p — n d i^d = ^) ^ s small, was achieved in the first r = 1 step of the present 
calculation. 

The above discussion shows in particular that the regime of n^ 1 ' small corresponds to a positive dnf jdx ps 
/5x at x small. In other words, the crossover for x small is absent at r = 2 in the parametric regime 
A pd 2 ~ 2t pd 2 ~ ApdKpl ^ 4 *pp- The condition dnf H) /dx ps with n ( d H) of Eq. (J22J) bears some resemblance with 
the corresponding MFSB calculation. As already mentioned, dn d MFSB ^ jdx = at fixed A p d/t pp , t p d/t pp was shown^ 
to occur at small ^ FSB for appreciable n d MFS (e.g£ for A p d ~ 4t p d S> \t pp \ when w 1/2 and n d MFSB ^ w 3/4). 
This indicates that dn d /dx w might also be achieved at the r = 2 level. However, expectedly, // ( ^ 2 - ) — becomes 
appreciable for close to 3/4. Such large variations of fi^ with r signify that the absolute convergence of the 

( 2) 

present theory becomes poor close to the crossover. The precise numerical analysis of the condition dn d 'jdx ps 
would thus be unreliable and one can use instead < 1/2, n d < 2/3 as a reasonable estimate of Xc7 from the 
x > x"cs side. 

(2) (2) 

It should be finally pointed out that D~ and given by Eqs. (|15ll6ll8p are not sensitive to the omission of 
the Cu-0 (anti)commutation rules. E^ of Eq. (TTJ]) is determined by B^' and F^ 1 -* which in turn are determined 
by Z?^ 1 ' and its / <-y b symmetric counterpart <fi^ . The latter are determined by , which is insensitive to the 

(2) (2) 

Cu-0 commutation and, therefore, so are Dl > and Tl '. This line of reasoning can be extended even further by 
' ' k k ° 

assuming that E( r ) ~ *iC , which is named^ 3 - the Ud = oo non crossing approximation (NCA). Since the NCA 

is / b symmetric it is A- independent and satisfies « 1 "on average" for rid small. It can therefore be called a 

"conserving" approximation^ According to previous numerical slave boson NCA^ and finite Ud DMFT calculations^ 

such corrections are expected to spread the d 10 od 9 disorder effects all over the coherent spectrum. However, these 

calculations^ start from the n d = 1 HF ground state of &-fermions which is not locally gauge invariant, in contrast 

to the n = ground state of Eq. ([T]) and those results should be therefore taken with some caution. Although 
Eqs. (|16ll8p show that the dynamical d 10 od 9 CT disorder induced by t p d is an inherent property of the large Ud 
Emery model in the regime A p d > 4|f pp | under consideration, further analysis is required to determine the full temporal 
(u! ± irf) structure of the higher order CT terms, before comparison with experiments. Most importantly, the present 
analysis indicates that the intrinsic dynamic d 10 od 9 charge disorder effect is "expandable" in terms of rid small, i.e. 
that it does not affect the Emery particle propagation strongly in the vicinity of the Fermi level, provided that rid is 
small. It should be kept in mind however that a strong electron-phonon coupling, omitted here, can slow down^ 2 - the 
d 10 <R-d 9 CT disorder and bring its effects to the Fermi surface. 



VII. VERTEX CORRECTIONS TO SINGLE PARTICLE p- AND rf-PROPAGATORS 

The lowest order diagram for E that is not taken into account by the NCA is shown in Fig. [5] As in previous 
low order diagrams, the arrows of time are carefully included to show the flows of energy throughout the diagram. 
Fig. [5] gives a local vertex correction to E^ 2 ) related to the local effective b — f interaction. It is clear that when the 
anticrossing of p-bands and the d-state is important the p-propagators in Fig. [5] have to be interpreted as the pdp 

— 2 f 1) 

hybridized propagators t pd T-, in the strictly iterative scheme. Such diagrams take into account that two permanent 
p-holes cannot hop simultaneously to a given Cu-site, one by the creation of the fo-boson and the other by the 
annihilation of the spinless fermion because the Cu-site carries either boson or a spinless fermion at any one time. 
Those contributions are allowed in the NCA, and thus have to be subtracted from Eq. (fH?)) . However, when rid is 
small the number of unwanted coincidences is expectedly much smaller than the total number of incoherent processes 
and this correction will therefore be ignored here. 

Further on, the two "squares" Ag in Fig. [7] which involve b and / propagators cannot be reduced to the self energy 
renormalizations of one of the propag ators or which appear in the diagrams of Fig. [7] The latter thus 

also extends beyond the NCA. E^ 3 ' depends on k and so it contains a nonlocal component. In this respect, we 
emphasize the appearance in Fig. [7] of two particle-particle or particle- hole pdp-\mes connecting two i^Ag's. In the 
SFT, the single particle lines can carry arbitrary spins. The corresponding elementary bubbles open singular particle- 
hole and the particle-particle correlation channels. This phenomenon is the seed of the pseudogap effects of the 
magnetic (SDW), charge/bond (CDW/BOW), and Cooper pairing correlations within Eg . Those terms and their 
extension to higher orders are therefore of utmost physical importance, in addition to the nonmagnetic r = 1,2,3 
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Figure 6: (Color online) Skeleton contribution (a) to T,^ generated by effective two-particle b-f interaction (b). 



A „ A 



b) 



A 



Figure 7: (Color online) Skeleton diagrams which contribute to irreducible El : (a) through the particle-particle channel; (b) 
through the particle-hole channel. 



terms describing in particular the local, dynamic d 10 <H>d 9 CT disorder. Apparently, the explicit nonlocality of I]! 3 ' 1 is 
useful for comparison with the NCA and the DMFT resultsM"-^ 



Figure 8: (Color online) Skeleton diagrams for two &«-»•/ symmetric effective p-particle interactions. 

The terms for El 3 ' additional to those shown in Fig.[7]are easily found on noting that the and Bf^ propagators 
must occur / <-> b symmetrically in El , in order to preserve the local gauge invariance. This corresponds to the 
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replacement of Ag in Fig. [7J by Ag + Ag (independent of A) for each spin, shown in Fig. [BJ The four- leg vertices shown 
in Fig. [5] describe the effective interactions between the pdp hybridized particles. When the fc-dependences of the 

— 2 (1) 

triangular vertices are associated with the pdp propagating lines, e.g. with t pd , the effective interaction between 
the pdp propagators is local (i^Ao), where Ag' b are the convolutions of two Bf ^ and two propagators, e.g. 



,4 A 6 _ ,4 2e rf - CJi - UJ 2 

L pd ly — L pdri I o- \ ' ^ ) 



Here, the ui s denote the external frequencies which run counterclockwise around Ag in Fig. [51 with u)\ in the upper 
left corner, so that cji + oj-2 = W3 + W4. Remarkably, ip d Ag, independent of A, is the product of four poles. Fig. [5] show 

that t£ d A ' are kinematic interactions, in the sense that one particle has to wait for the other to leave a given Cu 
site before it hops to this site. The causal nature of this result is associated with the — Hr\ position of the poles in 
Eq. (23). 

Analogous kinematic terms were first invoked 4 ^ in the T = (single band) theory of transition metals with sizeable 
Ud and recently extended^ to high T c cuprates in the single band limit. Eq. (|23|) thus represents the generalization^ 
of this concept to the multiband Emery model in the small A pc j, = 0, Ud = 00 limit. It should be noted in this 
respect that, in the regime n d < 1/2, the effective interaction t pd A is reasonably small itself. Indeed, upon assuming 
it to be small, one affects only the hybridized states on the Fermi level. Consistently, the important values of t pd Ao, 
small to moderate, are of the order of 



t pd A ~-^- <tpd . (24) 

^dp 

The reason is that, in the A 2 pd > 2t 2 pd > A pd \t pp \ > At 2 pp regime A dp = e d - p {1) > t pd for p,^ < uj vH (e d - uj vH 
depends only on t p d and A. p d and, e.g., is linear in t pd for small A p d). The SFT theory therefore converges quickly 
towards approximate local gauge invariance even when the processes of Fig. [7] are taken into account. 



Figure 9: (Color online) An infinite order expression for Eg(w); the first local term corresponds to the NCA. 



Actually, the diagrams shown in Fig. [Jj combined with Fig. [8l are the only skeleton diagrams for Eg based only on 
the vertex t pd Ao. Assuming that the effect of t pd Ao on is dominant (e.g. enhanced by nestin g 26 i 41 i 48 ) it is of some 

interest to construct the corresponding r — > 00 A-independent expression for 5]!°°^ on neglecting the off-diagonal b — f 
interactions of Fig. [S] The full infinite order partial sum, independent of A, is formulated in terms of the renormalized 
value of Ao, A and as follows. First, one introduces the infinite order slave-particle Dyson propagators in Ao 
instead of the bare ones, denoting the result by A . In Fig. [9] this is depicted by the open square. Second, as suggested 
by Fig. [7J and analogously to the usual single band perturbation theory, one joins Ao with the fully renormalized 
four-leg vertex A, shaded square, using the pdp-bnbble formed from two t~ d T^'s, denoted as arrowless in the hgure. 

Convoluting the result with t pd ^j. one finally obtains E^(w), which is depicted symbolically in Fig. [9] It is noteworthy 

that the nonlocal nature of Eg then corresponds to the nonlocal nature of t~ d Tg. The backward arrow in the second 
term of Fig. [S] is shown in order to stress that, in generalizing the skeleton diagram of Fig. [7J there must be a pdp-lme 
which starts and finishes in Fig. [9] with receding rl ^ of Eq. ((SJ, so involving the occupied fc's. The second term in 

Fig. M thus vanishes for 1 + x = while the first reduces to E^ ) = D^°>, as the exact Eg should. On the other hand, 
for 1 + x = the vertex t pd Ao of Fig. [5] remains finite, in contrast to that of Fig. [SJd. The approximation of Fig. [5] 
is thus nearly exact in the low density (electron doping) limit. This will be taken up again in Sec. IIXI while, in the 
following, we continue to discuss the approximation of Fig. [9] with the rid small in mind. 
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VIII. SPIN DISORDER AND RVB FLUCTUATIONS 



Figure [9] for omits the b — f interaction of Fig. [6] but not necessarily the "diagonal" b — b superexchange 

interactions stemming from Fig. [TU] which can be partially included in Fig. via the renormalization of single 6-particle 
propagators. Those effects, as well as the effects beyond Fig. [HI are conveniently discussed using the spin- flip correlation 
function x^((f, t) = — {i/N){T^2t bt b? , 6g f (t)bt _ (t)) also related to the spin flip dpd — dpd correlation function 

Xdd(q, t). Indeed, it is long know n 59 ! 67 from the slave boson theories that the super exchange interactions appear as 
interactions between two auxiliary spin carriers b a . The reason is that the spin-flip operator on the Cu-site = c|c^ 

maps on blbi in the slave particle theories. It means that two holes of opposite spins on an intermediate O-site 
can hop simultaneously, each to one of two empty Cu neighbors. Notably, the effective interaction of Fig. |TD] is 
purely antiferromagnetic if t~ d T^'' < with t pp — is used, otherwise the effective interactions between next-to-next- 
Cu neighbors are generated. The internal square in Fig. [10] is thus, quite schematically, of the order of the squared 
number of occupied pdp states, (1+x— rid) 2 . Although such an interaction is appreciable for |ar| small, the zeroth-order 

particle-hole correlation xib amon g the auxiliary spin carriers vanishes i.e. Xbb with B^ = B^ ^ of Eq. (|5|). The 
external lines in Fig. [TOl should be therefore associated with f?W = B 1 * 1 ^ + B^ K of Eq. (|A.3[) . The corresponding 
elementary correlation function, proportional to B^ * B^\ 



- 1 V" 1 2 b g 

<-bb V w ) — AT 7. /rwA (1) 

N ^ oj -oji L ){k) + e db > H 

i ^ 5 r^j /n^^ -f/~i (- 1 i i~t lfln/linrr r^nA i £i t i cirri lrfiKlo fr\i- m * ^ ^ 
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is therefore linear in a; n d to the leading order, i.e. negligible for n d small, compared to the corresponding 
correlations Xpp between the pdp-propagators, discussed in Sec. [X] Moreover, in sharp contrast with xpp i Xbb ls 
dispersionless and incoherent in the range p^ — > oj > ojm — s d \ , i.e., Xbbi 1 ^) l& rea ^ around oj = 0. Notably, 
the spin-flip Xbb describes the local dynamic spin disorder on the Cu-sites and is the only one renormalized further 
by effective interactions. 



A' 







Figure 10: (Color online) Skeleton diagram for b — b interaction, which generates the superexchange. 



Indeed, the analog of the genuine superexchange interaction J p d(q,oj) (i.e. RKKY J, J', J") of the permanent 
the Cu-sites appears (together with mixed terms) when the internal spinless fermion lines in 
Fig. EH are taken as F (1) = F^< + F^> (/•<->& symmetrically with Eq. IpPjl ) rather than as F^ = F^< : the 

empty pdp states are then involved too. This interaction itself is thus schematically given in powers of n d , i.e for 

n^ 1 ' small it is small with respect to the effective interaction tp d Ag between the pdp-propagators (not surprisingly, just 

the contrary was suggested^ 7 - upon coming from the x < x cs , n d = 1 , Mott side) . However again, in sharp contrast 

to tp d Ao, the superexchange interaction involves the .F , W > disorder, / <h> b symetrically with Xbb ■ ^ n the small n c i 
limit the overall effects of the magnetic interactions between the 6-particle-hole propagators come thus out as small 
and incoherent with respect to those of the interaction tp d Ao between the pdp particle-hole propagators. Notably, 

X bb (q,oj) resulting from the convolutions Xw, (<*>) * JpdiSi^) * Xhb ( w ) can then be interpreted as describing a dilute 
gas of RVB singlet o 2 ' 76 within the uncorrelated spin disorder of Eq. . 

(2) 

It is appropriate to add here some qualitative considerations about what can be expected when n d approaches 
1/2 for x > x cs which makes J p d and t^ d Ao comparable. This corresponds first to the evolution from the gas to the 

liquid of RVB singlets. Next we note that incoherencies of xll ( w ) ana - Jpd(q>u) are related to the dispersion of the 
L-band. In higher order calculations required for x w x cs the local AF RVB are likely to progressively localize the 
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pdp-states within the C11O2 unit cell, i.e. it becomes better to take empty pdp states as nearly dispersionless. In 
particular, intercell t pp should then be practically renormalized out,— which increases rid with respect to . The 
resulting Mott-AF coherence of the RVB singlets amounts to their short range AF ordering. This in turn plays in 

favor of increasing the number of empty pdp states i.e. tends to increase the spin polarized rid further. On balancing 

(2) 

such Mott-AF effects against the steady decrease of n d with decreasing x towards x cs from above one thus expects 
so to reach the saturation dn^/dx = at small x Rj x cs . Such an expectation is corroborated by the x < x cs slave 
boson calculation^ which starts perturbatively from the = 1 Neel state and shows that rid in the x = Mott-AF 
phase is in low orders larger than rid hi the small x > x cs metallic phase, all parameters of the Emery model except 
small x being kept equal. 

(2) 

Although this kind of analysis may be extended in principle to evaluate the magnetic correction to x cs we limit 
ourselves here to the observation that the low order SFT is consistent for x > x cs even when rid is somewhat larger 
than n d ■ 

IX. MODIFIED SLAVE FERMION THEORY (MSFT) 

The present approach derives the effective pdp-pdp interaction t^ d Ao in the small to moderate A p d regime, where 
it represents the dominant effective interaction. Its T = causal properties are explicitly related to local gauge 
invariance. Its shortcoming is, however, that it does not take into account the Cu-0 anticommutations. This is 
reflected in the spin structure of the effective interaction tp d Ao of Figs. Hk,,b. Due to spin conservation in the triangular 

vertices of Fig. [5J the effective interaction t^Ao of the hybridized t~ d T^ particles is the same in the singlet as in the 
triplet channel. Thus the assistance of slave fermion particles is required to forbid the double occupation of the Cu-site 
in both spin channels. However the triplet scattering in the original Hamiltonian is expected to vanish identically, 
irrespective of the value of Ud {Ud — 00 included) due to Cu-0 anticommutations. It is in this aspect that the theory 
which accounts for the Cu-0 anticommutations a posteriori, requires a modification of the SFT. 

We start by emphasizing that adding intermittently the first p- or d-particle (1 + x = 0) has to lead to the 
coherent single-particle hybridization irrespective of the value of Ud and the anticommutation rules, as was discussed 
in connection with Fig. |H1 The SFT reproduces this result by reducing Sr to So 5 i- e - Tr = and Dr = dQ\ 
i.e. only the coherent p — d hybridization is retained (irrespective of whether the anticrossing is important or not). 
Likewise, one finds B = BM and F = F<-°\ in a manner analogous to the phonon propagator in the standard polaron 
problems i21 

Next we turn to the N(l + x) = 1 "bipolaron" case. This is the fundamental problem in which the two requirements 
of no double occupancy and the d — p anticommutation rules come into play simultaneously. There are two ways 
to approach the problem of two particles in the system. In one approach, a single N(l + x) = 1 fermion, say with 
spin f, is placed in the k = [7r,7r]/a state of the L-band and the SFT is applied to the propagation of the additional 
p (or d) particle with spin either j" or J,. This is covered in the SFT by the slight adaptation to the fact that the 
spin in the ground state is now unpaired (time reversal symmetry is broken). The alternative approach involves the 
simultaneous creation, on the top of the empty p-band, of two particles in the singlet or in the triplet state. Obviously 
both approaches must coincide in giving the essence of the "bipolaron" physics. 

Let us thus start by considering the cohabitation of two p-particles in the SFT as a scattering problem. The fact 
that singlet and triplet scatterings of Fig. [5] are equal in the SFT bears some resemblance to the text-book scattering^ 
of two distinguishable particles of spin 1/2 by a local interaction. We see, therefore, that by treating the /-, b- and 
p-particles as distinguishable, ignoring consequently the d-p anticommutation rules, the SFT forbids double occupancy 
for both triplet and singlet scattering equally and dependently on band parameters Ed and t p d- 

On the other hand, with the original Emery Ud Hamiltonian it is the Pauli principle which forbids two particles in 
the triplet configuration to occur simultaneously on the Cu-site. The triplet scattering by the local interaction and its 
effects therefore vanish identically. By contrast, the scattering of two particles with opposite spins by Ud is finite and 
in the Ud = 00 limit has to correspond to the effect of forbidding the simultaneous double occupancy of the Cu site. 
In the conventional strong Ud coupling multiband perturbation theory only this effect can give rise to the effective 
singlet repulsion of two p-fermions independent of Ud but depending on the band parameters, t p d in particular. 

Returning to the slave-particle formulations this means that, in the MSFT, we have to remove by hand all dia- 
grams in which the two internal Mines in vertices of Fig. [5] carry the same spin projection and keep only those with 
two opposite projections. tp d Ao of the diagrams in Fig. [5] then represents the n d = counterpart of the n^ ' = 1 
resul t 59 i 63 i 66 i 67 for the superexchange J p d — in the Ud = 00 limit for A p d large. In the standard scattering 

problem^ with a local interaction, the a posteriori antisymmetrization of the wave function of two spin 1/2 particles 
makes the triplet scattering amplitude vanish, simultaneously doubling the unsymmetrized singlet scattering ampli- 
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tude. In the present case, Fig. [5J describes, in second quantization, the singlet scattering of two p-particles, which are 
not completely distinguishable because anticommuting between O-sites is taken into account. The problem therefore 
requires additional care. 



a) 




b) 
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Figure 11: (Color online) Irreducible b f+ / symmetric renormalizations for ' quadratic in t p d involving: (a) the local 
renormalization of the 6-propagator and (b) the local renormalization of the /-propagator. 

In determining the normalization factors of the scattering terms in the MSFT it is important to note that the 
SFT/MSFT self energy corrections, which amount to treating one particle in the "bipolaron" problem as permanent, 
are obtained by closing one of the two p-lines of Figs. [5Ji,b. Actually, rather than dwelling upon the "bipolaron" case, 
where the MSFT is to be exact, we can immediately address the limit of small. This reflects the fact that the 
necessary and sufficient condition for the fast convergence of the present theory is that is small, which is of course 
fulfilled, especially for 1 + x small (in this latter case for arbitrary band parameters). The resulting AE^ 1 ) is shown 
in Fig. [IT] Clearly it corresponds to the leading contributions to E^ of Eq. (|T4| given by 

AE (la » + A£( 1{, )< = —(AB[ 1] * F< + B> * AF™) , (26) 

where AB ( X 1} = (B^) 2 ^ (uj) and Af| 1} = (F^) 2 ^ (w) are shown in Figs. [3J and UJ and the spin is conserved in 
each triangular vertex. Fig. [TTJ shows that the effective interaction Ag' b between the hybridized particles occurs even 
in the leading contribution AS]' 1 ) to E^. In Fig. [TTJ Aq' 6 are, however, dissolved in the bubble renormalizations of 
and F^~\ This is unrelated to the d-p anticommutation and we must therefore introduce the requirement that 
A£W = A£$. In the SFT the spin factors of the first and the second term of Fig. [TTJ associated with /3^\u>) 

and 4> x (w), are equal to 1 and 2 respectively. The time and the Pauli structure of the theory multiplies the first 
term by — 1. The overall result is the subtraction — 1 + 2 = 1 of the two relevant terms in Eq. (|26p . In the MSFT 
the contribution of Fig. Ilia is removed by hand because it involves two bosons of equal spin, while the second term 
(Fig. Illb ) carries the spin weight 1 because the contribution with two bosons of equal spin is removed from the sum 

over the spins. This amounts to the (-F^ ) 2 <^ /2 renormalization of the F^ line. The overall result is again 1, i.e. 
AE (1) = A£$, provided that the prefactor of A b in Fig. [5b, i.e. the singlet scattering in the MSFT t^ d (A% + Aft) 
is itself taken equal to the singlet scattering from the SFT. As required, the d-p (anti)commutation rules are then 
entirely irrelevant in A£W and AEv, . 

Next, the relation A£^ = AE^ 1 ' is to be extended to the Dyson theory which removes the double poles (F^) 2 
and (B^) 2 from the above discussion. As previously discussed in the context of the SFT, I]' 1 ) is insensitive to the 
p-d (anti)commutation rules. The relevant question is therefore whether or not the choice = consistent 
with A£W = AE^ } , is unique, since the structure -B\ * (F^) 2 <p^ /2 of AE^ may suggest otherwise, despite 

the fact that AE$ = AE (1 >. However, expressing in terms of F^ ' and B^ ^ (/ •O- b symmetrically) allows the 
MSFT to benefit from the rapid convergence of the SFT towards the local gauge invariance and, in particular, to the 
LSR. In other words, by taking E^ = E' 1 ) the singlet and triplet Aq' 6 are taken again as dissolved, as in Fig. [TTJ 

in the bubble renormalization of B^' and F^ . At the r = 1,2 stages of the Dyson perturbation MSFT it is thus 
irrelevant whether the double occupation of the Cu-site with equal spins is forbidden by the slave particles or by the 
Pauli principle. 
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As already argued within the SFT, this line of reasoning can be extended within the NCA to T,^ r = E^ X ', i.e. 

(r) 

to Tj, , because here the interaction vertices are also absorbed in Dyson bubble renormalizations of all the /-, b- and 
p-propagators involved. By definition, the NCA omits all except bubble renormalizations, and is therefore the same 
for the SFT and the MSFT. The first contribution in Fig. M thus obviously corresponds to the "exact" SFT/MSFT 
NCA for the p-, /-, and ^-particles. 

It is instructive to examine the modifications related to the lowest order SFT expression beyond the NCA, namely 
to E' 3 ) discussed in connecton with Fig. [7J for r = 4. The introduction of the d-p anticommutation rules by the 
suppression of the triplet scattering of the physical fermions replaces (A-independent) E^ 3 ) by (A-independent) E^f ■ 
The difference E^ 3 ) — E^- is necessarily non-zero, because the interactions ^(Ag + Ag) in Fig.[5]cannot be dissolved in 

bubble renormalizations of the slave particles. E^ 3 ' — s|? has a non-local and a local component. While the non-local 
component in question is related to the breakdown of the Cu-0 anticommutation rules and should be discarded, the 
local component is involved in the SFT n d ^ + = n d ^_ anticommutation rule, as understood in Eq. (J9j) . The latter 
is thus affected in the MSFT beyond the r = 4 level of iteration even though the Dyson structure of Eq. ^ is 
maintained. On the other hand, the MSFT is taken to keep the empty d 8 state of Eq. ^ at "infinite" energy Ed + Ud 
as available for the creation of an additional fermion on the Cu site. This can be extended to r > 4 in Eqs. $Q and 
([7]) for E^ and rj^, obtained by replacing E( r_1 ) by E^ ^ and omitting the contributions of the triplet scatterings. 
In order to satisfy the equality n d r ^~ = n d + for r > 4 one should thus allow for the additional spectral density at 

the energy e d + U d on the t < side. This means that in Eq. © o- { [f 1]> differs from a^- l > of the SFT for r > 4. 
In other words, it can be reasonably expected that the MSFT is locally gauge invariant on average and, to a good 

(r) 

approximation, Q MR ~ 1 if the / -<-> b symmetry is maintained in the calculation of the /- and b- propagators. One 
can even expect that the MSFT is a better "conserving" approximation^ than the NCA because it has the advantage 
that it generates incommensurate magnetic correlations even for finite t pp . The asymptotic local gauge invariance is 
however conserved only if the omission of the triplet scatterings between any two p<ip-particles, corrected for the lack 
of the local spectral weight by Eq. (|9]), renders the ground state exactly antisymmetric. The rigorous proof of this 
assertion is lacking so far. Nevertheless, the MSFT has some advantages over the DMFT based on small cluster o 54 i 55 
of Cu02 units, because this latter is appropriate only for the description of commensurate structures. Unfortunately, 
the DMFT approach which allows perturbatively for incommensurate magnetic correlations^ is at present restricted 
to the single band Hubbard model. 



X. LADDER APPROXIMATION FOR COHERENT MAGNETIC CORRELATIONS 

The present theory based on the Ud — > oo Hamiltonian of Eq. (H} generates the d 10 od 9 disorder in the r = 2 
single particle propagators and the singlet pdp particle-particle and particle-hole correlations for r = 4, as shown in 
Fig. [71 Since the effective interaction t pd Ao turns out to be repulsive for long times its primary effect concerns the 
spin flip electron-hole correlations in addition to the d 10 -f->-d 9 disorder. Although the d 10 od 9 disorder occurs in the 
single particle propagation in lower order than the spin flip processes, the leading r — 2 disorder contribution falls 
far from the Fermi level and the subsequent r = 3 disorder term expands in terms of Ud, assumed small according 
to Eq. (|2~4")) . It should be noted in this respect that the choice of the iterative renormalization scheme is not unique, 
i.e. that other appropriate prescriptions can be devised according to the value of x and the band parameters under 
consideration. E.g., the NCA emphasizes the d 10 f>d 9 charge-transfer disorder effects and neglects the magnetic 
correlations. In contrast, the choice of the free and in Fig.[7j3 together with the hybridized pdp propagators 
t pd T\ ' obviously favors the strong coherent SDW correlations for x w x v h- Notably, this flexibility in the selection 
of the perturbation subseries has its physical counterpart because magnetic coherence is expected to compete with 
the d 10 <H-d 9 disorder, and win when it is sufficiently strong. This suggests the investigation of Fig. [JJd using the free 
i?(o) ) B (o) and t^T^ propagators. 

In this context, it is interesting to examine the pdp-pdp magnetic spin-flip susceptibility xsdw i n the coherent 
limit. The first relevant spin- flip correlation is associated with the free pdp-pdp bubble \pp ( < ?! CJ )- For lu small it is 
dominated by the intraband contribution Xpp(Q>ui)- The latter should be distinguished from the dpd-dpd bubble with 
the intraband contribution Xddi^ 00 )- Both functions exhibi t 26 i 48 the same singular behavior in the reciprocal space 
when pf^ falls close enough to the band van Hove singularity but differ in the associated spectral densities. 

Xpp is further renormalized, in particular according to Fig. fTS] where i 4 d Ag appears explicitly as the effective 
singlet interaction between the r = 1 pdp propagators^ An analogous contribution comes from ^Aq. The actual 
calculations are simplified by the separability of the interactions ft d A^' b , like in Eq. (|2"3")l . The eight poles involved 
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Figure 12: (Color online) tp d Ag renormalization of the magnetic spin flip susceptibility; spin (vertical arrows) is conserved in 
triangular vertices. 

are split in 4 + 4 sets. The real part of the integrals is approximately proportional to (x PP ) 2 tpd^o hr agreement with 
the estimate of Eq. (|24[) (note that this results in cancellations of numerous factors between the squared spectral 
weights z^ of Eq. (fTU)) and tp d /A dfl ). Apparently, the retarded nature of i 4 d Ao does not play an essential role in 

this result, (x pp ) 2 tpd-^o thus increases the bare Xpp - This latter can be large at finite x due to nestin g 26 i 48 and the 
ladder approximation then leads to strong low energy SDW correlations xsdw, 

Xpp 

Xsdw ~ Tjy j (27) 

1 - UdvXpp 

where ~ 4i 4 d /A 3 ^ [ of Eq (f2~4")l . Equation ([2T)l corresponds to the simplest vertex renormalization (which omits to 
single out the Umklapp component 4 ^ of Ud^) in the symbolic equality shown in Fig. |pj 

Such a distribution of the roles of the propagators and interactions contrasts with the small Ud theory for xsdw, 
where the interaction Ud involves the dpd bubbles Xdd ■ So does the ansat a 37 ' 79 Ud — > J p d(cosk x + cosk y ) with 
Jpd = 4ip^/Ap d that is traditionally used in the approach from the MFSB side. Note the essential role of the negative 

sign associated with the cos term for k = [tt, n]/a and questioned in Ref. IzE Equation ([27]) resolves thus the long 
lastin g 37 ' 76 ' 79 controversy of what replaces J p d = 4i 4 d /A 3 d of the Ud = oo Mott-AF limit in the Ud = oo metallic 
limit. The answer is that the quantity in question is the repulsion Udp ~ 4tp d /A.^ between the pdp propagators. 

In cuprates this replacement occurs progressively by doping through x cs as discussed in Sec lFXI All contributions 
to the magnetic susceptibility xsdw, those stemming from xlt of Eq. (|25|) and the cross terms should in principle be 
treated on equal footing with xsdw near above x cs . However, xsdw dominates the coherent fluctuations. As long 
as \(j, — u} v h\ is large, so that the "nesting" effects in Xpp (<fj w ) are weak, the spectral weig hts zi L) of Eq. CEOJ may be 
important in this respect. Under the additional assumption that \fi — uj v h\ is small xsdw of Eq. (|27p becomes large 
at uj small and its singular behavior determines then alone the low frequency behavior of the system. 

A similar line of reasoning, invoking the effective, instantaneous repulsion between (simply renormalized) single 
particle propagators, with additional corrections, can be applied to the discussion of many observed properties of 
cuprates^ ' 38 ' 39 ' 41 It is therefore of some interest to discuss briefly how the present theory deals with some salient 
experimental results. 

XI. COMPARISON WITH EXPERIMENTS 
A. ARPES and NQR 

ARPES experiment o 3 ' 4 ' 30 ' 31 measure the dispersion properties of the exact electron propagators in the reciprocal 
space which are related by the sum rule to nd and n p determined from the NQR data£ They show unambiguously 
that for x > x cs fi falls in the broad band, below the Cu site energy id- As in the MFSB case^ the nonmagnetic 
(unreconstructed^) band structures with "Fermi arcs on large surface" are observed for x > x cs . They are first fitted 
here to determine two renormalized band parameters, conveniently t p d/t pp and A p d/t pp , while t pp is chosen finally 
to fix the overall energy scale. The observed Fermi energy fi is conveniently measured by lj v h — (J- where uj v h is the 
position of the vH singularity. In parallel, the NQR measurements determine the absolute values of nd{x) with a few 
percent accuracy due to the ambiguity in the evaluation^ of the Madelung contribution to the electric field gradients. 

The measured band parameters are identified here with those calculated in Eqs. (fl"6|) . (fT5|) and (fl~7|) for r = 2, 
namely with t^~ d /t pp and A^/t pp (keeping in mind that e p and t pp are not renormalized). For comparison with 
experiments the latter should include iteratively or selfconsistentlyr 9 - primarily through the site energies e^ 1 ' and 
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e p , the variation of the Madelung potential due to the variation from n d to n d . Those parameters, together with 

— M^ 2 ^; determine rS 2 ^ and x e ff of Eq. (|22p . In principle the bare parameters t p d/tpp, A p d/t pp and jJ-^ can be 

found from t^J /t pp , A p x J /t pp and x e jf under assumption that renormalizations are weak. In practice this procedure 

implies a heavy numerical work required by the inversion of Eqs. (| 1 5|) and (fllj)) . Thus, instead of calculating rff from 
ARPES the latter is estimated from NQR measurements of rid and n p . In this spirit we discuss below two different 
r = 2 regimes, associated respectively with lanthanates and YBCO. 

The example of an excellent, essentially two-parameter fit obtained from the three band model is given in Fig. [5] 
for x = 0.07 LSCO. Typical values of the band parameters satisfy A 2 d > 2t pd as A p d\t pp \ > it pp , the regime with 

relatively flat Fermi arcs. The observed /x corresponds to n d 2H ^ m 1/2. Upon doping x = x v h ~ 1/5 (1/8 in 3 - 1 - LBCO) 

/i traverses^^ the vH singularity. The Luttinger sum rule is qualitatively satisfied^ - indicating that n^ 2 /2 as n d 2lnc ^ 
cancellations in Eq. (|2"2")l lead to x e // ~ x. Indeed, good band fits can then be obtained in function of x, as already 

noted before^ 9 - in the MFSBT. The renormalizations t p X ] /t pp and A p X J/t pp of band parameters by are significiant, 

but NQR gives^ rid < 3/4 and a positive drid/dx w 1/4. Although consistent with the decrease of t pd /\t pp \ and increase 

( 1) (2) (2H) 

of A p J /\t pp \ upon doping, we feel that the resulting difference between rid ~ n d and n d somewhat stretches the 
low order theory. 

For YBC0 6+l5 NQR gives 5 . n d ~ 1/2, dn d /dS « 1/10 and x « 8/5 for x > x cs w 0.05, while the observed 8 ^ 
Fermi level ii falls well below the vH singularities (of the band doublet) at uj v h — A* ~ 0-25 eV which does not 
depend much 8 - 2 ^ on S. The observed A p d/\t pp \ and t p d/\t pp \ are smaller than in lanthanates, towards the regime 
A pd « 2t pd « A p d|ipp| ~ 4i 2 p exhibiting pronounced arcs of the unreconstructed 8 -^ Fermi surface. Consequently, 
the r = 2 renormalizations of band parameters are weaker and n d is smaller than in lanthanates, consistently with 
large u>^ — The negative incoherent corrections to x e ff of Eq. (f2"2"j) tend then to be large compared to n^ 2 /2. 
Notably, x e f / ps x — n d sa means that doped holes go at average mostly to the highly degenerate set of disordered 

d 9 states^ which is consistent with d(uj[^ - (i&)/dx 0. The ARPES result - fi)/d& « 0, traditionally 91 

named breakdown of the Luttinger sum rule (assuming a reasonable a; ~ 6), can thus be brought into broad agreement 
with the theory. 

While ARPES provides d(fj, — L0 v H)/d5 for x > x cs , it is generally believed^ 3 - that XPS in cuprates gives the absolute 
shift dfi/d8 itself. It turns out^ 3 . so that in metallic YBCOg+s (5 « 5x), d/j,/dS < (in the hole language) and tends 
to vanish when optimal doping 5 ss 1 is approached. On the other hand, for metallic LSCO, the shift d^/dS (8 ~ x) 
is found 8 - 3 - small all along. This issue is discussed here by setting conveniently 

= - + & - \ ^Af 2 + m p f (28) 

where e^ 1 ) is the average site energy = (e d + e p )/2. As easily seen such expression cannot reconcile the ARPES 
data and the XPS analysis unless the variations with x of the in plane Madelung field produced by the negative 
charge x of dopants are taken into account. In particular, the component of the Madelung field which is uniform 
within, and among, the Cu02 unit cells produces the variation of with x, while A pd {A pp = e px — e py = 0) 
are less affected. The estimation of the average Madelung field is simple when the negative charge is distributed in 
homogeneous layers of surface density x/a 2 . The negative shift S M e of (conjugated^ 3 - to rid + 2n p ) of the order 
od xVpdc/a for lanthanates, unscreened within the plane, tends to cancel out the positive variation of (i due to the 

filling of the Cu02 plane with x holes at fixed This agrees broadly with the XPS analysis^ 

Notably further, dfj,/dx oc (6 (rid + 2n p )) -2 <C 1 indicates 5 ^ a propensity of the normal metallic phase to inhomo- 
geneous phase separation in the plane. Macroscopic charge redistribution is however hindere d 43 ! 44 ' 59 by long range 
Coulomb screening. It is thus appropriate to consider the charge transfers over short distances, among or within the 
Cu02 unit cells, as will be discussed next. 

B. Magnetic orderings and lattice deformation structures 

As mentioned in the Introduction, various experimental methods, 1 ^ including NMR, 3 ^ 3 - STM^ X-ray2i and 
quasi-elastic neutron scatterin g 18 ! 20 in particular, reveal entangled ; 24 ' 25 (quasi-)static magnetic and lattice deformation 
structures, which break, or restore as a checkerboard, the D4 symmetry of the Cu02 lattice. They are characterized 
by the Bragg spots which occur at qsDW and <fo respectively and, assuming the satellites to those leading harmonics, 
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are usually called stripes i 21 ' 23 At small < x < x cr w 0.05 go and qsDW lie on the zone diagonals but for doping 
x > x cr they are rotated by tt/4 and become collinear with the main zone axesJ^ The relation 

go + 2q SDW = G (29) 

is satisfied all alongi 14 i 21 i 23 In addition to stripes, LBCO-like cuprates exhibi t 19 i 27 i 28 the commensurate LTO/LTT 
lattice phase transition. 

The present approach attributes diagonal stripes to the weakly doped Mott-AF phase. Indeed, it is likely that 
a coherent t,t', J, J' models can roughly 7 ^ explain such stripes, though we are not aware of any previous to ours 
attempt^ 8 - to relate t', J' etc. to parameters of the Emery model. On the other hand the present theory of the metallic 
phase should produce the collinear stripes or checkerboard configurations which satisfy q$f u + 2q c g^ w = G. 

/ 
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Figure 13: (Color online) Umklapp contribution due to two SDW's at qsDW ~ G/2, resulting in momentum transfer qo. The 
vertex denoted by the dot depends on the symmetry of the intracell fluctuation or phonon which carries the momentum go 
away. 

When dealing with this issue we realize first that external frequencies in the involved correlation functions can be 
taken to vanish (adiabatic approximation), i.e. it is then possible to construct the corresponding Landau-like theory. 

In contrast to X^(<f, 0) of Eq. (|2"5|) Xpp\qsdw ~ G/2,0) of Eq. (f2"T|) exhibits a logdike singularity at lo pz which 
is the strongest when \lo v h — ^\ is small enough, as it certainly is in lanthanates, as we choose from now on for 
brevity. Importantly, Xpp (<fi w ) ma y then presen t 48 i 51 i 79 peaks for incommensurate gl^y in the vicinity of G/2 even 
when t pp is finite. This can be taken to generate the Landau theory in which the q^°/j w SDW is the primary order 
parameter. The lattice deformations are driven by two such SDW's through the Landau invariant shown in Fig |13l 
obtained on assuming that all external frequencies are small. It can be easily seen by direct inspection of Fig. [TBI that 
when \uj v h — mI i s small the main contribution comes when three pdp propagators which form the triangle fall into 
the vicinity of the vH points. Insertion of the Bi g -like vertex on the top of the triangle in Fig. [13] is then unessential, 
i.e. the two produce a large Bi 9 -like O x /O y CT with the wave vector q^ oU which satisfies Eq. (|2"9"1) . 

while the Logdike CT is small. The induced Ai 9 -like CDW on the total charge of the Cu02 unit cell is also small for- 
go smal l 43 1 44 and moreover screened^ 4 - by long range Coulomb forces. Unscreened, the latter is sizeable for arbitrary 
<fo and, in particular, interacts^ 3 - by the e-p interactions with LTT modes. Moreover, the LTT tilts are entangled, 
either via O x /O y CT or directly by ionic interactions^ with the q^ oU shears of the Cu02 plane. A sizeable cfe° u is 
necessarily associated with linear coupling of the lattice displacements, such as is the oxygen LO phonon^ to the 
mixture of O x /O y , Cu/02 CT's and CDW (including the bond density waves not discerned here). Let us note on the 
semi-quantitative side that even when q^ oU is close to the Cu02-lattice tetramerization 27r[l/4, 0]/a, the omission of 
Coulomb screening in the straightforward MFSB calculation^ of stripes overestimates the q^ oll -CDW with respect to 
the q% oll -O x /O y and Cu/0 2 CT's. 

The observed q^ oU is large with respect to the inverse SDW correlation length ; 14 ' 18 ] 20 i.e., the formation of stripes 
does not interfere to the lowest order with the behavior of Xpp^((f, u>) at q « 0. Kexpp' L \q, 0) is logarithmically 
large at (J, ~ uj v h (% ~ x v h) in lanthanate a 3 i 30 not only at qsDW but also at q = 0. This corresponds to the 
Jahn- Teller (JT) splitting of the vH singularitie a 1 ' 43 ' 44 rather then to their nesting, relevant for the SDW case with 
\^>vH — Ml small. Such JT splitting produces the homogeneous q w intracell ^>\ g -O x /O y CT n px — n py , while 
the i?2 S -contributions and the overall Ai 9 -CDW are absent by symmetr y 43 ' 44 Clearly^ the dynamic H>i g -O x /O y 

CT is accompanied by local though not circular— currents. \pp (<? = 0i w ) (without spin flip) has no direct vertex 
renormalizations by large Ud (in particular the Coulomb interaction V pp between two neighboring oxygens, may be 

invoked for further enhancement of Xpp )■ The q = O x /O y CT is coupled by the quadratic e-p coupling to two LTT- 
tilts with Q = 27r[l/2,0]/a each, since the linear e-p coupling then vanishes^ 3 - by symmetry. Such LTT dimerization 
(which produces the splitting A pp = e px — e py of the oxygen site energies, conjugated^ 3 - to n px — n py ) plays the role 
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of the primary order parameter— of the first order LTO-LTT phase transition. The latter gives rise in particular to 
the new Bragg spot at 2Q = 27r[l,0]/a, as it should ) 19 ' 29 accompanied by the stripe leading harmonic at q^ oU . 

The above description agrees remarkably well with the observe d 24 ' 25 SDW- LTO-LTT structure of the striped phase 
in lanthanates, including the first order LTO-LTT phase transition in LBCO. It emphasizes the importance of sizeable 
O x /Oy CT n px — Tipy between two oxygens and is thus consistent with appreciable charge sharing between Cu and 20 
invoked here. In contrast, it is obvious that those results and in particular the simple O x /O y JT effect, can hardly be 
explained by the rigid t — J model which keeps the average (oxygen) occupation of the upper band low, proportional 
to x. 



C. Inelastic spin-polarized neutron scattering 

The inelastic spin polarized scattering of cold to thermal neutrons is usually taken to measure the exact w), 
associated with the flip of the overall spin in the Cu02 unit cells. This fixes the relative phases of xVki Xpd an< ^ Xsdw 
in x (QiW). The large scattering for oj around and above 0.15 eV occur s 18 ' 20 for u>{q) which differs clearly^ from 
the AF-magnon parent dispersion 6 ^ with J p d ~ 0.15 eV. In addition the dispersion is considerably broadened and 
presents the absolute maximum of intensity for q^ g DW (uj) on zone diagonals. As concerns Xsdw(q, w) of Eq. (|27[) . its 
properties result from the interplay in Xpp \Q} aJ ) °f "small" energy scales, ui, lo v h~H, t pp and t p d(q—G/2). A systematic 
analysis of Xsdw(& w ) is currently under way^ but early numerical calculations^ - of x' (<T, T) for reasonable values 
of band parameters indicate that Xsdw{q, w) may be enhanced at appropriate values of (fsDW^ )- The corresponding 
energy transfers uj are also possibly affected by xlt associated with Eq. (|25|) and its RVB extensions. Further on, 
serious departures from the Heisenberg AF-behavior occur at lower frequencies, especially below u)hg ~ 50 meV where 
the so called hourglass dispersion is observed^ This indicates that Xsdw(q, w), which includes the single particle 
(pseudo)gapping,— due here to the SDW itself, is important in the frequency range ujhg > ui > u p h where ui p h are the 
bare frequencies of phonons involved in static stripes. In this range xsdw{q, w) decouples from the lattice (but not 
from the O x /O y CT in Fig. I13[) and an interesting open question is whether it by itself tends to break or not the D4 
symmetry of its ladder approximant xsDwiQ,^)- The previous analytic^ and numeri c - 51 ' 79 calculations for such ui 
show that xsd w(Q> cj ) presents peaks at q^° D l w on the main axes, in agreement with the observed 7r/4 rotation of the 
absolute maxima in the neutron scattering. Finally, for ui < ui p h ~ 12 meV one observe s 18 ' 87 the conic dispersion of 
spin-lattice density waves, expected in generalization of the present approach which takes the SDW at tfg DW {u) ~ 0) 
as the primary order parameter of the stripe formation. 

D. Magnetic pseudogap 

Beyond the unreconstructed Fermi surface, such as that in Fig. [5j the ARPES measurement s 4 ' 31 ' 88 reveal in lan- 
thanates the pseudogap around the antinodal (vH) points of the Brillouin zone, while in YBCO-like (and electron 
doped) materials the strong spectral density variations^ 9 - are observed along the unreconstructed Fermi surface, to- 
gether with shadow bandsi 80 ' 81 These structures are usually associated with the magnetic energy scale of about 
50 meV, reminiscent of u>h g associated above with the hourglass dispersion. Such ARPES structures were discussed 
befor e 42 ' 60 ' 90 in some detail, assuming that the propagating hole emits and annihilates a damped soft AF-paramagnon 
described by a phenomenological x^(<f f=s G/2,u) which couples to the propagating hole by an adjustable coupling 
constant. Although incommensurability and couplings to the lattice are omitted in x , those approaches reproduce 
observations remarkably well. 42 ' 60 ' 90 In particular the adjustable coupling constant of those calculations is associated 
here by Fig. [7] and Eq. ([TT])) . with the microscopically derived Ud^ of Eq. (f2"T|) involved also in Xsdw- This lends a 
semi-qualitative support to "paramagnon" approache s 42 ' 60 ' 90 to the magnetic pseudo gap. 

E. Raman, optic and soft-X-ray measurements 

The compelling evidence for the crossover picture proposed here comes from the Raman and optic experiments. 
We notice in this respect on the expected magnetic scale up to 0.5 eV that the energy of the sharp two-magnon 
resonance observed^ 6 - at x — in the Bi g Raman data for lanthanates is about 0.35 eV. This Raman activity arises 
from the coupling^ 7 - of the two-magnon resonance of the Heisenberg model on Cu sites with super exchange J p d (via 
O-sites) to the q = Bi g quadrupolar O x /O y CT within the Cu02 unit cell (and further to photons). Measurements 
show clearl y 16 ' 17 that the two-magnon resonance shifts downwards and fades progressively away for x > x cs as 0.05. 
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This is necessarily related to the energy shift, broadening and loss of the spectral density of magnons involved in the 
two-magnon resonance. The continuity through x cs suggests that the broad structure observed for x > x cs reflects 
the short range AF-RVB fluctuations generated by Eq. (|25|) . Indeed, the contribution of two SDW's each described by 
Xsdw(q-, w) coupled according to Fig.Q2]to the quadrupolar q = qo O x /O y charge CT within the C11O2 unit cell is not 
expected to produce large Raman response at the energies as high as 0.35 eV. Actually the pdp triangle in Fig. [T3"lis 
dominated for \uj v h — n\ small by three propagations in the vicinity of the vH points, i.e. this contribution is expected 
to be large only for lower frequencies, and then at finite qo(u) 7^ 0. Rather than with Fig. [IS] the Raman q = Bi g 
response at frequencies below Uhg ~ 50 meV is thus dominated by the q = O x /O y CT associated to the Jahn- Teller 

XpJ (0, w) with pseudogap corrections. So far however there is no direct evidence for the coherent quadrupolar exciton 
(in the whole Raman range of frequencies) and the O x /O y CT fluctuations are observed through their coupling to 
phonons. In YBCO this concern s 16 ! 91 e.g. the q = B lfl phonon at 340 cnr 1 and in lanthanates the Q LTT-phonons, 
coupled in pairs to Raman photons via the intracell quadrupolar O x /O y CT-fluctuations. This assignation agrees 
qualitatively with sharp and two phonon-broadened features observed respectively in x ~ 0.1 YBCO and lanthanates 
by low frequency Raman 9 - 1 - scattering. 

Similarly, the optic conductivity of conducting lanthanates exhibits a hump^ 3 . at w w 0.5 eV for < x < 0.05 
which shifts downwards, broadens and eventually fades away upon doping through x cs ~ 0.05. This hump may well 
correspond to the transitions of x holes involved in the Zhang-Rice singlets between the bonding and nonbonding 
states on oxygens^ 3 - It is reasonable to think for x > x cs that the (downshifted and broadened) 0.5 eV optical hump 1 - 3 - 
is the RVB remnant of the Mott-AF phase, in parallel with that observed^ 6 - around 0.35 eV in the Bi g Raman 
scattering, while the rest can be explaine d 44 i 45 by the interband optic transitions, such as those of Fig. [5j of the 
weakly renormalized three-band model. Indeed, no resonant optical structure related to the 3-band model appears 
in the 0.5 eV - 50 meV range of our r = 2 calculation. Notably, the crossover from the Mott-AF state to the 
metallic state is accompanied in such approach by the shift of chemical potential with x from that associated with 
small x Zhang-Rice structure to that appropriate for the moderately renormalized 3-band model. As concerns the 
conductivity at frequencies below uj hg w 50 meV, it is described here basically in terms of (no spin flip) xg?(0,w) 
with nonmagnetic vertex corrections, includin g 44 ! 45 long range Coulomb screening. This picture may also take into 
account the self-energy renormalizations of the pdp-propagators associated with magnetic (pseudo 42 > 60 i 90 ) gap^ 4 - scale 
and perturbative vertex corrections associated to nonmagnetic Umklapp scatterings, such as that arising from the first 
and second neighbor Cu-0 interactions V p d and V pp . The vH antinodes are however unimportant for low frequency 
conductivity, because the carrier mobility vanishes there with or without a pseudogap. In contrast to that the vicinity 
of nodal points is not^ affected by the pseudogap at least for sizeable t pp , giving rise to a conduction picture, which, 
in the simplest possible low T approximation, is associated with the Fermi liquid Umklapp T 2 law in the dc resistivity. 
Interestingly, by using sum rules appropriate for coherent electron-hole excitations^ the described picture for infrared 
conductivity and Raman activity can be brought in good agreement with the Hall effect measurement s 13 ! 33 ! 92 in weak 
magnetic fields, while the corresponding quantum oscillations^ are currently under consideration from the present 
point of view. 

The soft-X-ray absorption is associated with the creation of electron-hole pairs, with the hole created in the deep 
core state. Actually, an early observation^ 5 - of the crossover at small x cs was accomplished using this method. The 
interpretation^ 5 - of the data was carried out using the Zhang-Rice t — J limit^ 3 - for all x of interest. For x = and 
x > holes are annihilated respectively in the lower level and in the upper level (narrow ZR band) , using the hole 
language. While we keep the same picture for x < x cs , we associate here the X-ray absorption in the x > x cs metallic 
phase with the transition of the hole from the renormalized lower L-band to the core level. The difference AE between 
the (low order) hole energies for the x < x cs and x > x cs regimes is ; 49 ' 61 i 68 in particular, according to Eq. Q28p. 



taking for simplicity that x cs — + and t pp = 0. This takes into account that, unlike in the metallic phase, the intercell 
hybridization is forbidden in the x — Neel phase. Appreciable renormalization of the Cu-0 hopping in Eq. (|17p . 
rather than of the average site energy s and the CT gap A p d, consistent with Fig. O makes SE negative. The latter 
is also approximately equal to the negative shift in the luminescence frequency, as originally observed^ 5 - Although 
such a trend is satisfactory it is appropriate to mention that, when the variations of na — 2n p (conjugated^ 3 - to A p d) 
between the x < x cs and x > x cs regimes are sizeable, ea and e p also vary due to the concomittant variation 5 -^ 9 - of 
the Madelung field. The variation of A^ in the Neel (N) phase with respect to the value of for x = + is 
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while e is essentially unaffected for x cs = + . The Madelung constant olm is dominatec^ 9 - by the first and second 
neighbor interactions V p d and V pp , w 2.14. With repulsive 3 ^ V p d ~ 2 eV and Jn^ rj 0.1 — 0.2 the difference of 
excitation energies is thus increased further already in the iterative (rather than selfconsistenl^) Madelung scheme. 
Those observations reconcile, for x > x cs , the ARPES data where the hole is created at the Fermi level and the 
measurement of the soft-X-ray absorption^ where the hole is annihilated at the Fermi level. In both cases, the Fermi 
level is consistently assumed to lie in the renormalized but still broad L-band of Fig. [5] 

F. Relation to superconductivity 

In our Ud — > oo limit the superconductivity is itself a coherent, translationally and locally gauge invariant* (double 
occupancy of Cu forbidden) state. It is therefore also expected to smoothen out the intrinsic disorders to some 
extent, but apparently it competes with the SDW and O x /O y CT coherences. Extraordinarily, the superconductivity 
disappears entirely in the commensurate LTT phase of LBCO for x w x v h ~ 1/8 (/i ~ w„g) when the q = O x /O y 
CT gets frozen by its coupling to the commensurate 7r[l, 0]/a ~LTT/(s X x — £yy) lattice deformation accompanie d 19 ' 29 by 
its q^ oU satellites. This provides a striking evidence for the importance in high-T c superconductivity of D4 symmetry 
(approximate in YBCO7) involving the in plane oxygens. Actually, the superconductivity may even be enhanced 
when the Bi g O x /O v CT corresponds to the dynamic D4 symmetric (intra- or inter-band) quadrupolar exciton since 
the latter can contribute to the bosonic glue^i^" 3 - However, as already mentioned, a coherent exciton is not observed 
within the frequency range accessible to B lg Raman scattering^ In contrast, the 340 cm -1 B lg phonon in YBCO 
coupled to superconducting order is well know n 16 ' 91 for a long time. 43 ' 93 Notably in this respect, the superconductivity 
is stronger in YBCO than in the lanthanum cuprates, where the zone boundary 7r[l,0]/a LTT/(e X2; — e yy ) phonons 
are evidently much heavier. Such issues apparently require further in-depth investigations, including the possibility 
of superconducting (rather than AF) space coherence among the RVB singlets.— At present we can only infer directly 
from Fig. [5] that the effective repulsion Ud^ = 4£ 4 ld Ao of Eq. (|23p may generate the Coulomb pseudo potential 11* in 
the pdp-pdp Cooper ladder channel^ which is not prohibitive to high-T c superconductivity, especially in materials 
like YBCO, where rid is smaller— than in lanthanum cuprates. 

XII. SUMMARY 

Careful consideration of the overall experimental situation indicates the existence of a crossover between the insu- 
lating and metallic state for x — x cs (A p d,t p d,t pp ), which is quite small in actual cuprates since the model parameters 
obey A 2 d > 2t 2 d > A p d\t pp \ > 4i 2 p . The x < x cs regime can be associated with the Mott-AF or short range AF 
regime, while the x > x cs regime corresponds to the modified band picture with magnetic correlations. Such a dis- 
tinction explains the persistence of magnetic correlations through the doping x cs and contrasts with the description 

T A 



SF 


/ SB 




SB 


SF 


fSF] 


SF 




/ AF 




AF 


stripes 1 








/ + extrinsic disorder 






CHARGE + SPIN 


\IO 


1 CHARGE 










DISORDER 1 


LTO 


| DISORDER 












LTT 




-1 1 


4 1 1 — 







1 - 


— 1 — 


J 1 x 



— I I — 

YBC0 7 LBCO 



Figure 14: (Color online) Cuprate behaviors associated to slave particle theories; SF stands for MSFT and SB for the anologous 
slave boson perturbation theory; shown are stripes, commensurate LTO/LTT instability and the spin plus charge disorder of 
the d 10 «d 9 type. 
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of the hole doped metallic regime by the rigid t — J model for all x>0 of interest. It is self-evident that the modified 
band regime also applies to the metallic state for x < (the MSFT is exact for x = — 1), where the enhanced stability 
of the insulating AF phase is related^ to the pinning of doped electrons to the dopant sites (to the extrinsic disorder). 
Once the x > x cs time-dependent perturbation theory is shown to posses the local gauge invariance asymptotically, 
it is employed in the usual way, namely by summing up selectively chosen subseries, which emphasize the physically 
important features. The question left open here, whether or not the omission of the triplet pdp-pdp scattering is 
sufficient to antisymmetrize the SFT a posteriori, thus becomes of limited practical importance. 

This is summarized in Fig. 1141 which also shows the main physical features observed in cuprates, at least briefly 
mentioned here, namely the intrinsic CT d 10 -H-d 9 disorder, the stripes and the commensurate LTO/LTT instability 
driven by the coherent, purely intracell O x /O y CT. In the metallic x > x cs regime all these are well described by 
the present theory, together with many other properties. Remarkably, the optimal doping in highest-T c cuprates, 
represented in Fig. [T?] by YBCO7, falls in the range x cs < x < x v h, while in lanthanates it is close to x = x v h where, 
moreover, superconductivity is suppressed^ by the commensurate LTO/LTT lattice instability. 

In conclusion, the MSFT contains features which qualify it for further development in the context of the high T c 
superconductivity, just as do the corresponding moderate to strong coupling multiband theories in terms of finite Ud, 
some of which were discussed herein. 
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Appendix: Evaluation of and E (1) 

The purpose of this Appendix is to calculate explicitly B^\ and E^ 1 ) ~ B^ * Fp defined in Eq. (fM)) . B^ 
is given by 

= - 1 ^TT— = B IP(» - A), (A.l) 

uj - e d - A + irj - pi (io) 

where /3P'(u>) is defined by Eq. (fTTj) . Here +irj keeps the memory of the pole in B^ while the poles in /3P\uj) are 
on the opposite side of the w-axis. In the usual procedure, one uses the well known equality 

P- ± mS(x), (A.2) 



X =P IT] X 



and the N — > 00 limit to calculate (3^ ' (uj) in integral form. It is then important to retain +irj in Eq. (|A.1[) because 

Eqs. (fTTj) and (| A.2|) show that Im f}^ (uj) vanishes strictly everywhere except on a finite segment of the w-axis. Im B^ 
turns out to be negative on the whole Rew-axis as appropriate for bosons. The information about the side of the 
w-plane taken by the poles is, however, lost in this way. Analogously, ImF^ changes sign twice which is far from 
the conventional (Fermi liquid) behavior of the fermion propagator. Additional care in the evaluation of , 
id, Aj[ p in Eqs. (fT5")) and (fTljj) is thus required. 

To this end we consider Eq. (|A.1[) as an equation involving a discrete set of poles, however dense. According to 
Eq. (fTTj) the propagator B^' contains the group of N/2(l + x) (receding) poles in the upper half plane and one 
(advancing) pole in the lower half plane. This means that B[l) can be written in the form 

A oj - E d b - A + irj A 



B 



(Ah)/ w(L) 

vE ; A),* , , — (A-3) 



(1)< = 

u-u,l L \k)-\ + 2ir, 
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where the unique advancing pole B^ ' is separated from B x ' . Its position idb + A is shifted from Ed + A according 
to £< db = £ d + P\~ (w = £< db + ^) ! anc ^' expanding /3^\uj) around + A it is supplied with the spectral function 
1 + nY/(o;)/2 written in the way convenient for the boson on the t > side. We keep A in order to check how it 
evolves through various steps of the actual calculations. is represented as the sum of the receding poles shifted 

to the positions oj^ L ' = uj^ L \k) + AwJ L '(fc) with the appropriate spectral densities r y^ b \oj). 

It is, however, out of reach to evaluate analytically the energy shifts and the spectral functions for each of the 
1 + N(l + x)/2 poles in Eq. (|A.3j) for finite N and therefore we turn to the N — ► oo limit of Eq. (|A.3j) . One starts by 
using Eq. (|X2|) in Eq. (JTTJ) for ^ (w) 



It is worthy of note that Eq. (|A.4[) exhibits the unshifted positions and the spectral weights of the N(l + x)/2 poles. 
Those poles are grouped in the energy interval [w^u*] where ujf n = ljm + A and uj^ = + A are the energies of 
the lowest and the highest occupied pole in the L-band shifted by A. uim is the shorthand notation for the k = [w, 7r] 
energy in the L-band and fjS 1 ^ is the HF chemical potential. 
We assert next that in the N — > oo limit 



1 ImP^\uj) 

t uj - e d - A - (RepM(w))* + (Im/^H) 2 



= Im^ 1)< (w) (A.5) 

using the expression (|A.4[) for Imf3 x \uj). To show this we compare Eq. (|A.5|) with Eq. (|A.3[) . All sums for 

Im/3^~\— w)'s in Eq. (|A.5[) should not be converted into integrals simultaneously. Neither should they all be kept 
as sums since the ^-singularities should not appear simultaneously in the numerator and denominator of Eq. (|A.5p . 

Indeed, according to Eqs. (|A.3|) and (|A.2j) Im B 1 ^^ can be expressed as the sum of N(l + x) /2 (5-functions, associated 

with the values of uj [L \k) shifted to Cj^ L \k) = uj( L \k) + Auj { b L \k), 

lmB^<(w) = ^T^M^^ - 4 L \k) - A) . (A.6) 

k 

The shift &.uj\^\k) is of the order of 1/N in the N —t oo limit but, nevertheless, it has to be taken into account 
when defining this limit in Eq. (|A.5j) . This is most simply achieved by performing the fc-integration for Re (3^ and 
in the denominator of Eq. (|A.5j) while still keeping Im/3x in the numerator as a sum over N(l + x)/2 poles 
at the positions infinitesimally shifted to 6j( L '(k). Such a procedure eliminates "on average" the J-functions from the 
denominator and leaves the infinitesimally shifted ^-functions under the sum in the numerator. 

The comparison of Eq. (|A.6|) with such a modified Eq. (|A.5|) allows us to identify 7^ Ab ' (w) of Eq. (|A.7[) as 

t 2 z {L) 

k oj - e d - A - (Re 0™ (w)) 2 + (Im (w)) 2 

for each of the N(l + x)/2 poles after letting Aw£ (A) — > 0. Whenever associated with 5(w — lu~) in its numerator 
^^ b \uj) can be replaced by 7jj ^(^r) by setting w = uih = + A in the " A-averaged" denominator of Eq. ()A.5[) . 
It is obvious that such an ImB x vanishes strictly outside the interval [w^u;^] due to the fact that the 1/AT shifts 
Aw[ L '(fc) are neglected after the N —> oo limiting procedure. Importantly, is needed in Eq. (|A.7[) in addition 

to Jm/3 X to take care of the "average" spectral density of the "unshifted" receding poles. Inserting "f^ b \uj) of 

Eq. (|A.7|) into Eq. (|A.3[) now determines . 
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Obviously, nw(u) = e ( ib + A) in B 1 ^^ is nothing but n^}, independent of A. The structure of n^ 1 ' is however more 
transparent when calculated on the t < side from n^) , 

= -2 / ImB^^duj . (A.8) 

Eq. (| A. 8[) can be evaluated in two ways. One can insert the result (|A.7[) for ^^(cu) in Eq. (|A.8|) . carry out the 
w-integration and then the fc-integration. Alternatively one can transform the numerator of Eq. (|A.5I) immediately in 
the integral over k and then carry out the w-integration in Eq. (|A.8|) . The result is the same, i.e. the two integrations 
in question are interchangeable in the N — > oo limit. Noting that Ari[^ of Eq. (|12p corresponds to omitting /3r (w) 
in the denominator of j{ Xb \co), it is immediately apparent that n£ ~ ^ n b = n d^i provided that is small. 

The analogous b O / symmetric procedure can be carried out for the spinless fermion F^\uj) = F^^uj — A), as 
already mentioned in connection with Eq. (fT5)) . In associated with AB^ 1 ' and AF^ 1 ' of Figs. [3] and [5] one sums 

over two bosons according to Eq. (|A.8[) , and AQ^' = 0. However, differs from unity due to the asymmetry of the 
roles of the spin factor 2 in and F^\ In particular the receding pole in F^ shifts from A to A — 2/3^ (u> = e^+A). 
Similarly, 7 1 A /) (a;) differs from 7 1 A b) (w) of Eq. (|A7f|) by the appearance of the factor 2 in the denominator. It is then 
clear that ps 1 for small when n£ ps Au^' and n^ 1 ' pa A^ . 

Finally, using By~ of Eq. (|A.3|) and its 6 O f symmetric counterpart F^ one finds that 

4 1) =£d + 3^ 1) (w = 4 1) +A) (A.9) 



eP appears in E' 1 ^ of Eq. (fT5|) which corresponds to the convolution of only two fro/ symmetric poles B 1 ^^ and 

Fi . S^ 1 ^ is taken here in the N o oo limit with the short hand notations (2h = + A and 0)2 = u> — ef +Sd + A. 

As easily seen S( 1 )< is, like E^ 1 ^, independent of A. Equation (|A. 10|) provides the closed f -> 00 expression for 
Ap p, in Eq. (fTrJ)) . It can be easily turned into integration by using Eq. (|A.2[) . 
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